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Chapter 1 


Ha tea*® & ® m. 


This document describes the International Reference Ionosphere 1990 
(IRI-90). It is intended as a guide and handbook for the experienced as well 
as the new IRI user. IRI models are established by a Joint working group of 
the Committee on Space Research (COSPAR) and the International Union of 
Radio Science (URSI). The composition of the COSPAR/URSI Working 
Group on IRI is listed on the inside front cover. 

Close ties exist to the URSI Working Group G.3 on Ionospheric Modeling (C. 
M. Rush, U.S.A., Chairman) and G.4 on Ionospheric Informatics (B. W. 
Reinisch, U.S.A., Chairman) and to the SCOSTEP (Scientific Committee on 
Solar-Terrestrial Physics) Working Group on an Aeronomical Reference 
Ionosphere (R. W. Schunk, U.S.A., Chairman). 

IRI describes monthly averages of the electron density, electron 
temperature, ion temperature, and ion composition in the altitude range 
from 50 km to 1000 km for magnetically quiet conditions in the non- auroral 
ionosphere. 

Almost a decade has passed since the last comprehensive IRI handbook was 
published: IRI-79 was described in Report UAG-82 of the World Data Center 
A for Solar-Terrestrial Physics (Rawer et al., 1981). Meanwhile, the IRI 
model has been significantly improved with ground and space data collected 
in the seventies. Work and studies that led to several IRI updates were 
presented and discussed at the annual IRI workshops and are published in a 
series of issues of Advances in Space Research: Volume 2, No. 10, 1982; 
Volume 4, No. 1, 1984; Volume 5, No. 7 and No. 10, 1985; Volume 7, No. 6, 
1987; Volume 8, No. 4, 1988; Volume 10, No. 8 and No. 11, 1990. The list 
of contents of these reports is reproduced in Chapter 5. 

The IRI-90 handbook summarizes the most important improvements and 
new developments. It includes a text part, a portion that explains the IRI 
formulas and expressions, and finally a collection of figures generated with 
IRI-90. 
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An interactive driver program for IRI-90 was developed at the National 
Space Science Data Center (NSSDC) and is distributed (together with the 
IRI-90 FORTRAN computer code and coefficients) on magnetic tape, on 
diskette, and on line on the Space Physics Analysis Network (SPAN) and 
connected computer networks. The diskette version can be executed on 
IBM compatible personal computers (XT, AT), no math co-processor 
required*. 

Since 1988, IRI has also been part of NSSDC's Online Data and Information 
Service (NODIS), which can be accessed from a remote SPAN node: (1) SET 
HOST NSSDCA, (2) USERNAME=NODIS, (3) choose the GEOPHYSICAL 
MODELS option, and (4) follow the prompts and menus. You can (i) run the 
program, (ii) read the documentation, or (iii) get instructions on how to 
copy the IRI code and coefficients to your account. 

IRI-90 on tape or diskette can be ordered from — 

• For U.S. researchers: 

National Space Science Data Center 
Code 933.4 

Goddard Space Flight Center 
Greenbelt, Maryland 2077 1 
Telephone: (301) 286-6695 

Telex: 89675 NASCOM GBLT 
TWX: 7108289716 
FAX: (301) 286-4952 
SPAN: NSSDCA: -.REQUEST 

• For researchers residing outside the U.S.: 

World Data Center A for Rockets and Satellites 
Code 930.2 

(Address same as above.) 

Publication of this document was made possible through the support and 
encouragement of Joseph H. King and through the dedicated and untiring 
efforts of Miranda Knowles, Lynda Williams, and Maria Walters. 


* With co-processor, the execution speed increases by a factor of ten. Recompilation is necessary if the 
user has a system with co-processor and wants to take advantage of the faster speed. 


2 



Chapter 2 




Contents 


2.1 History of the International Reference Ionosphere 5 

(K. Rawer, L. Bossy) 

2.2 Electron Content Measurements and IRI 1 1 

(R. Leitinger) 

2.3 Morphology of Electron Temperature Anisotropy in the 

F-Region 23 

(K. -I. Oyama) 

2.4 Model Descriptions for the Ion Transition Heights 33 

(D. Bilitza, I. Kutiev) 

2.5 Models for Horizontal E- and F-Region Drifts 41 

(E. Kazimirovsky, E. Zhovty, M. Chernigovskaya) 




3 





2. 1 History of the International Reference Ionosphere 

K. Rawer 

Herrenstr. 43, 7806 March 

Federal Republic of Germany 

L. Bossy 

University Catolique de Louvain 

1 74, avenue W. Churchill, 1180 Bruxelles 

Belgium 

Radio-observations of the very first satellites yielded some very important 
results concerning the upper atmosphere. Orbital period determinations via 
radio-location showed that temperature and density variations in the upper 
atmosphere are much greater than expected. Therefore, the newly- founded 
Committee on Space Research (COSPAR) decided that a set of empirically 
based tables describing these findings should be established. This task was 
considered of such high importance that a commission of specialists was 
established to cany it out. This group in 1961 presented its findings under 
the name of COSPAR International Reference Atmosphere (CIRA). The 
collection was widely used by the scientific community. Taking account of 
the increasing amount of measured data followed up by an Increasing 
number of theoretical approaches, the continuous work of the CIRA Task 
Group culminated in the release of three more CIRAs In 1965, 1972, and 
1986. 

Realizing the success of CIRA, S. Bowhill proposed a few years later to 
establish a comparable reference for the Ionized constituents of the 
atmosphere, to be called International Reference Ionosphere. According to 
the original terms of reference, it should contain empirically based tables 
describing monthly median vertical profiles of the main parameters of the 
ionospheric plasma. Like CIRA, IRI should be usable as a standard reference 
for the design of experiments, the estimation of environmental and other 
effects, for checking theories, etc. The profiles should be provided for 
suitably chosen locations, hours, seasons, and levels of solar activity. In 
contrast to theoretical models, the IRI should combine approved 
experimental results so as to be a useful reference with no dependence on 
theoretical assumptions. 

To take care of this new task, COSPAR in 1968 established the "Task Group 
on IRI" and nominated K. Rawer for chairman. This group was given the 
task of promoting and coordinating the international efforts toward the goal 
described in the above terms of reference. In 1969, the International Union 
of Radio Science (URSI) decided to concur with COSPAR in the project. A 
preliminary set of tables (for two locations, noon and midnight, three 
seasons, and two levels of solar activity) was presented as an example to the 
URSI General Assembly 1975 in Lima, Peru. 
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The Task Group made two decisions at the very beginning of its work: 

( 1 ) Since most users have access to computer facilities, the results should 
be presented as a computer code rather than as a set of printed tables. 
Tables and figures could then be produced with this code. 

(2) The peak of the ionosphere should not be modeled Independently but 
should be determined by an existing computer code published in 1967 
by the Comite Consultatif International des Radiocommunications 
(CCIR). 

The CCIR code is based on a large set of measurements made by ionosonde 
technique at more than a hundred stations worldwide. This data base, 
however, contains large gaps over the oceans. In order to fill these gaps, 
Jones and Gallet (1965) had introduced a particular, coordinate-dependent, 
interpolation system that was bound to geographic coordinates. To give a 
more coherent picture, the CCIR version applies another, special latitudinal 
coordinate (MODIP) that was proposed by Rawer (1963). 

The terms of reference asked for four parameters to be represented in IRI, 
namely, electron density, electron and ion temperature, and (positive) ion 
densities. In order to reach consistency with the independently obtained 
electron densities, the Ion densities should be given as relative, not absolute, 
values. Consistency was also asked for between the plasma temperatures 
and the CIRA neutral temperature. 

As a general philosophy a data set, before being introduced into the system, 
had to be critically evaluated by the experts. Comparing results of different 
techniques was one way of checking data reliability, a powerful one as was 
found out. In fact, critical comparison of different measuring methods 
turned out to become a major IRI task. For example, a special meeting on 
"Methods of Measurements and Results of Lower Ionosphere Structures" was 
initiated by the group and was held in 1973 at Konstanz, Federal Republic of 
Germany. With the guidelines established during the meeting in mind and 
in broad international cooperation, a great effort was started to gather 
relevant data from different techniques. 

At the onset of IRI work, it had been expected that the amount of accessible 
and reliable data would be rather large for electron densities, much smaller 
for temperatures, and rather poor for ion composition. In general, this was 
found to be the case. However, even the electron density profile data 
showed rather important and unexpected gaps. In spite of the existence of 
an enormous amount of evaluated ionogram parameters, true height profiles 
were only available for a small number of stations providing by no means a 
worldwide coverage. The inversion technique, needed to obtain true height 
profiles, had been applied only at very few places at temperate latitude. 

Even there, some uncertainty remained because of the E-F-valley problem 
(Gulyaeva et al., 1990). Several stations run by U.S. institutions had 
produced so-called "composite profiles" obtained with an artificial "median 
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ionogram" (established month by month and hour by hour). Though the 
reliability of these results might be questionable, they were the only available 
input at lower latitudes. Also, the polar caps were de facto excluded because 
of the lack of specific information and because of the great variability 
encountered in this region. It was felt that under the extreme conditions in 
the polar caps the CCIR maps, which are based on monthly median station 
data, are unable to produce a representative picture. 

As for the topside, mainly two sources of information were available at the 
time, namely, incoherent scatter observations (unfortunately at few stations 
only) and topside sounder profiles. In situ measured plasma densities could 
be used only for checking purposes. A very large number of topside 
ionograms taken by the two Alouette satellites had been inverted to density 
profiles. Unfortunately, an inquiry showed that these were incomplete 
because the peak electron density value was regularly missing. Thus, 
regrettably, the many profiles in the archives were not usable for the IRI 
purpose. Later Bent and Llewellyn (1973) established an empirical model 
description based mainly on about 10,000 topside ionograms obtained by 
Alouette above North and South America. While the influences of solar 
activity and F-peak density are given by simple formulas In this model, the 
latitudinal variation is described discontinuously considering only three 
geographic ranges. S. K. Ramakrishnan produced a continuous description 
for the Llewellyn and Bent (1973) model. There remained, however, serious 
problems in the region around the geomagnetic equator where the original 
description does not admit an explicit dependence on latitude (Bilitza, 

1985a, 1986). 

For the height of the F2-peak, hmF2 , reliable data are obtained by the 
Incoherent scatter technique; however, only a rather small number of 
stations operate worldwide. Fortunately, in view of propagational 
applications, ionosonde stations had regularly determined a parameter, 
M(3000)F2, which has some relation with hmF2. A large and worldwide 
thesaurus for this parameter was readily available in the form of the CCIR 
numerical maps similar to the CCIR peak density maps. From a compilation 
of ionogram data (under some simplifying assumptions) Shimazaki (1955) 
had found a linear relationship between the peak height hmF2 and 
M(3000)F2, a relationship that was later improved by different authors 
(Bradley and Dudeney, 1973; Bilitza et al„ 1979). The Improved formulas 
take account of additional parameters which influence the hmF2-M(3000)F2 
relationship, in particular refraction in the E-region. With such relations 
the M(3000)F2 maps could be transformed into hmF2 maps. The results 
obtained by incoherent scatter had been taken into account for the 
improvement introduced by Bilitza et al. (1979). 

Particular problems were encountered in comparisons of measurements in 
the lowest ionosphere (below 100 km) that had been obtained with different 
techniques. The data were so widely different that a special symposium was 
held (see above) in order to reach an agreement about general guidelines 
(Rawer, 1974). It was stated that in situ measurements, when combined 
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with radio propagation measurements between ground and rocket, would be 
used as the primary input. A compilation of acceptable D-region rocket data 
was made by Mechtly and Bilitza (1974). For the nighttime lower E-region 
and valley two differing compilations were available from Maeda (1969, 1970, 
1971, 1972) and Soboleva (1972, 1973). Comparing both with Schumann 
resonances, H. G. Booker gave a strong vote in favor of the Maeda model. 

The full (vertical) profile of plasma density was described by a set of 
mathematical expressions, each valid in a certain height range. This rather 
complex system allowed a correct representation of the most important 
inputs, like the peak densities of the main layers. On the other hand, it was 
not well suited for "full wave" computations. For such computations (at 
extremely low frequencies), Booker (1977) proposed the use of fully analytic 
functions of a type that P. S. Epstein had defined in the thirties. His 
proposal was, in fact, used in the topside description but for good reasons 
could not be realized for the full profile at that time. 

At least by day the electron temperature can be much higher than that of the 
ions, the latter being also less variable. Ion temperature data observed with 
the retarding potential analyzer technique by Dumbs et al. (1979) were 
accepted as reliable inputs together with ground-based measurements made 
at a few places by the incoherent scatter technique. The in situ results were 
used to derive the latitudinal variation between the incoherent scatter 
stations. The difference to the lower CIRA neutral temperature values 
becomes quite appreciable above 400 km. Near the 100 km level agreement 
of both temperatures (as requested by COSPAR) was arranged for within the 
ion temperature formula. 

The same data sources were used for the electron temperatures. The 
latitudinal variation was mainly taken from a compilation of data measured 
aboard the AEROS satellites (Spenner and Plugge, 1979). Later Bilitza et al. 
(1985) improved the temperature model by introducing Langmuir probe 
data compiled by Brace and Theis (1981) from the AE-C and ISIS 1 and 2 
satellite missions. It is well documented that the electron temperature and 
density are closely anti-correlated in the daytime ionosphere (e.g., Bilitza, 

1975; Brace and Theis, 1978). Therefore, in the formulas used at first, the 
electron temperature was coupled with the actual electron density. Later it 
was found that for monthly medians a correlation between the temperature 
and density profile was not meaningful. Density and temperature profiles 
were, therefore, established independently from the relevant monthly 
median data. In case, however, actual density profiles for a specific time and 
location are at hand, it is recommended to apply a correction to the median 
temperature profile via an empirical (inverse) relation specified in IRI 
(Bilitza et al., 1985). 

The data base for (positive) ion composition was rather poor at the time. Most 
of the published rocket or satellite data had been given as absolute densities 
and did not provide the total plasma density needed to determine the relative 
ion composition. Further, a large set of simultaneous measurements with mass 
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spectrometers and a retarding potential analyzer (RPA) aboard the AEROS 
satellites had shown that while the spectral resolution is much better with 
mass spectrometers, the intensity indications are less reliable with the more 
Involved spectrometer systems than with the simpler retarding potential 
method. The RPA data could be reasonably well evaluated in three groups: (1) 
hydrogen and helium ions, (2) oxygen and nitrogen (atomic) ions, and (3) the 
group of molecular ions: 02 + , N 2 + , and NO + . Satellite RPA measurements 
(Dumbs et al., 1979) were used for the model at heights above 220 km. A few 
ion mass spectrometer (IMS) data (Taylor, 1972) were also considered in 
building the model. Between 100 km and 250 km the model relies on rocket 
IMS data compiled by Danilov and Semenov (1978), including American, 
Russian, and a few European measurements. 

The most difficult range to describe is below 100 km. Positive cluster Ions 
appear below about 92 km and are the most abundant species below about 
84 km. There were only three groups in the world that had developed 
techniques (with cooled spectrometers) as needed at heights at which 
collisions play an important role. It was not possible to establish a profile In 
agreement with all observations because cluster Ions are often destroyed in 
the spectrometer itself (by the electric fields that are applied in the 
measuring system). After some discussion, data obtained with a particularly 
designed instrument (Kopp et al., 1978: Kopp, 1984) were taken as 
reference. Only the total sum of all clusters was given since the number of 
individual species Is large and transitions between species is frequent. 

A similar situation was encountered with negative ions that occur only at 
altitudes below about 80 km. These are also clusters. Only one group had 
made measurements in this difficult height range (Arnold et al., 1971). 

Their average daytime profile was incorporated into IRI. 

After the first IRI was released in 1978 (Rawer et al., 1978a), it was 
critically tested with a wide variety of data. This testing period lasted about 
one decade. Since 1982, COSPAR in cooperation with URSI has organized 
yearly workshop meetings to discuss and improve the model (see Table 1). 

The computer code has been changed step by step and new features 
Introduced as they have become available. COSPAR has published the papers 
presented at these meetings in its periodical. Advances in Space Research. 

The proceedings of seven of these workshops have appeared in the years 
from 1982 through 1990 (see Table 1). More than 300 scientists are listed 
as authors of more than 230 individual papers filling a total of 1,458 pages. 

The tables of contents of these reports are listed in Chapter 5. 

Another result of these discussions is worth noting: For quite a number of 
unresolved questions, the existing data base was insufficient to give a well 
founded answer. As a reaction, special data analysis and collection efforts 
were undertaken in different countries. For example, the AEROS satellite 
team decided to organize the data reduction and evaluation scheme in such 
a way that a worldwide picture of electron and ion temperatures, and ion 
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Table 1: Workshop Meetings for Discussion of IRI 


Year 

Place 

Published 
Year Papers 

Publication 

Pages 

1971 

Seattle (USA) 

1972 

21 

Soace Res. XII 

1229-1335 

1973 

Konstanz (FRG) 

1974 

52 

See Rawer (1974) 

1 -460 

1974 

Sao Paulo (Brazil) 

1975 

5 

Space Res. XV 

295-334 

1982 

Ottawa (Canada) 

1982 

14 

Adv. in Space Res. 2. # 1 0 

181-257 

1983 


1984 

24 

Adv. in Space Res. 4, # 1 

1-169 

1984 

Graz (Austria) 

1985 

18 

Adv. in Space Res. 5, # 7 

1-112 

1985 

Louvain (Belgium) 

1985 

23 

Adv. in Sgace Res. 5, # 10 

1-130 

1986 

Toulouse (France) 

1987 

22 

Adv. in Space Res. 7, # 6 

1-127 

1987 

Nnvnorod (USSR) 

i 1988 

53 

Adv. in Space Res. 8, # 4 

1-251 

1988 


1990 

24 

Adv. in Space Res. 1 0, # 8 

1-132 

1989 

1 Abinadon (UK) 

1990 

17 

Adv. in Space Res. 10, #1 1 



composition could be obtained in a readily usable form. In the U.S.S.R. and 
in India several rocket campaigns were conducted to get a sound data base 
for an improved ion composition vs. height formula. 

Not all steps taken in the course of discussions were found to be favorable. 
While originally an interrelationship between electron density and 
temperature was implemented, it was found that, in spite of clear 
experimental evidence for individual cases, this is not justified when dealing 
with monthly medians. Another unfortunate decision (in response to a 
request by the competent URSI Commission) was made with the temporary 
introduction of Chiu's (1975) simplified peak description as a second choice 
besides the numerical CCIR maps. This action was withdrawn later when it 
became apparent that this (regionally helpful) description is not acceptable 
in many parts of the world. 

In their 1978 paper. Rawer et al. (1978b) noted a few open problems: 
checking by propagation experiments, comparison of elf propagation with 
full wave computations using IRI, comparison with ionospheric absorption 
measurements, extension toward greater heights, improved data base for ion 
composition, in particular at lower and at very great heights, and improved 
data base in the low latitude belt. Some of these problems were resolved; 
meanwhile, others are still pending. 
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2.2 Electron Content Measurements and IRI 


Reinhart Leitinger 

Institut fur Meteorologie und Geophysik 

University of Graz, 8010 Graz, Austria 

Abstract — The use of Ionospheric electron content data from Differential 
Doppler measurements for comparison with the International Reference 
Ionosphere (IRI) Is demonstrated. A data base of European electron content 
data exists in Graz and was used for the comparison. The data were 
obtained at several receiving stations from the signals of the Navy Navigation 
Satellites (NNSS). Their nearly polar and nearly circular orbit is at a 
nominal height of 1000 km. The data gained at Lindau/Harz, Federal 
Republic of Germany, were used, with calibration by means of latitudinal 
profiles of electron content observed from Graz, Austria, if available. The 
comparison was based on monthly medians of electron content for latitudes 
60°, 55°, 50°, and 45° north and local time intervals of two to three hours. 
Two levels of solar activity were considered: 0 £ R £ 40 (low solar activity — 
LSA) and 130 £ R ^ 170 (high solar activity — HSA). Apart from electron 
content, two other Integral parameters that can be derived from the IRI are 
discussed: (1) the equivalent slab thickness T and (2) an ionospheric shape 
factor rj. Both parameters are Important for the assessment of 
transionospheric propagation errors. 


2.2.1 Measurement Technique and Data 

Electron content and other integral parameters for the ionosphere from the 
ground to ceiling heights around 1000 km are Important both for 
geophysical and for engineering applications. Since actual measured values 
are not always available, one has to rely on models for many purposes. One 
of the areas of great practical importance Is the prediction and the 
assessment of transionospheric propagation errors as needed, for example, 
space geodesy, radio astronomy, time dissemination via satellites, satellite 
and space probe navigation. Only a few papers have been published that deal 
with the comparison of measured electron content with the IRI (e.g., 
McNamara and Wilkinson, 1983: McNamara, 1985; Bilitza et al., 1988). 
Almost all published comparisons were based on electron content gained by 
means of the Faraday effect on the signals of geostationary satellites. 

Modern measurement techniques for electron content are based on 
propagation effects observed on signals transmitted from satellites to ground 
stations. Several signal components are needed to separate plasma effects 
from geometric and neutral gas effects. The most important experimental 
methods are listed In Table 1 together with the signal components used. 
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Table 1: Measurement Techniques 


Technique 

Faradav Effect 

Differential Doppler Effect 


Alternate Name 

- 

Carrier Phase Difference 

Modulation Phase 
Difference 

Signal 

Components 

Left- and right-hand 
circular components of 
one carrier 

Two coherent carriers with 
large difference in 
frequencies 

Identical and coherent 
modulations on two 
carriers 

Frequency Bands 

VHF 

VHF/UHF, L-Band 

c:‘ VHF/UHF, L-Band 
m: t 0.1 .. . 10 MHz 

Satellites 

Beacon of geostationary 
communications satellites 

Navy Navigation Satellite 
System (NNSS), Global 
Positionina System (GPS) 

ATS 6 Radio Beacon 
GPS 

Carrier 
Frequencies 
f,. fs/MHz 

136 

NNSS: 150,400 
GPS: 1227, 1575 

ATS 6: 140,360 
GPS: 1227, 1575 

Reference 

Frequencv/MHz 

- 

NNSS: 50 
GPS: 1227 

- 

Sensitivity 

A^/ATEC" 

[ o /f10 15 m-*;i 

9.12 

NNSS: 9.24 
GPS: 15.5 

ATS 6: 2.11 
GPS: 0.13 


c - carriers. 

* m - modulation. 

M A* - phase difference; A TEC - electron content difference. 


Using geometric optics, the general expression for the received signal phase is 


4) = ^ n ds + 2nft 

T 

where f is the transmitted frequency, c the free space velocity of light, n. 
the refractive index, ds the ray path element, t the time, and the integral is 
taken from the phase center of the transmitting antenna at T to that of the 
receiving antenna at R. 

All three techniques record propagation effects in the form of phase 
differences: (1) the Faraday effect as the phase difference between the two 
circularly polarized components of a suitable satellite signal, (2) the 
Differential Doppler effect as the phase difference between two coherently 
transmitted signal carriers transformed to a common reference frequency 
by phase division, and (3) the Group Delay effect as the difference of 
modulation phase on two carriers. In each case, the phase difference 
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contains plasma effects only. The Faraday effect is approximately 
proportional to the integral over N where N Is the electron density and 
B l the component of the geomagnetic field strength along the ray path TR; 
the other two effects are approximately proportional to the Integral over N 
only. For signal frequencies greater than 100 MHz these approximations are 
very good and higher order effects can be neglected. In the case of the 
Faraday effect, the integral is evaluated assuming a constant value for the 
whole ray path. One usually takes the Bj, value from a fixed height. For 
geostationary satellites errors smaller than 5% are obtained with a fixed 
height of 420 km provided that the result is considered to be the electron 
content from the ground to a limiting height of 2000 km CHtheridge, 

1972). The contribution of the electrons above this limiting height is 
negligible because of the decrease of B L and of N with height. 

Accordingly, propagation effects on signals transmitted from satellites are In 
a first order approximation proportional to the electron content between 
satellite and receiver. This slant electron content can be transformed 
(approximately) into the vertical electron content at a path height h t , by 
multiplying It with the cosine of the angle between ray path and vertical 
direction (zenith) at h ( . It is common to take a h 4 of 350 km or 400 km. 
Model calculations have shown that h t = hmF2 + 50 km is a better choice 
with hmF2 being the height of the F-layer peak. Ionospheric electron 
content is often called total electron content (TEC) to distinguish it from 
sub-peak electron content or other partial contents. 

If the propagation effect gives a phase difference greater than 2k, the data 
are ambiguous and one has to use additional information to resolve this 2 rut- 
ambiguity. In many cases, the transmitted phase and/or the phase shifts in 
the receiving equipment sire not known, which means that one has to add a 
phase constant (an initial value) that can assume any value (not only 
multiples of 2k). Provided that transmitter and receiver are stable in 
respect to phase shifts, the initial value remains constant as long as there is 
temporal continuity in the received data. Therefore, the initial value 
problem is in general less severe for geostationary satellites than for orbiting 
satellites. In the latter case one has to find a new initial value for each pass 
of a satellite. For geostationary satellites the initial value acts as an additive 
constant because the projection from slant content to vertical content is 
done with a constant factor. During the pass of an orbiting satellite the 
projection factor changes and, therefore, the Initial value has a changing 
influence on electron content. For orbiting satellites only those electron 
content differences gained for satellite locations with equal zenith angle are 
free from the Initial value influence. For low and midlatitude stations and for 
satellites in nearly circular polar orbits, equal zenith angle is approximately 
equivalent to equal latitude differences to the latitude of the receiving 
station. 

Ionospheric electron content depends on time and on location. It is not 
possible to gain both dependences with one receiving station and with one 
satellite: The beacons on board geostationary satellites make it possible to 
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derive the time dependence of electron content along a fixed ray path (or 
vertical content for a fixed location). Satellites in nearly circular and nearly 
polar orbits in heights around 1000 km (e.g., the Navy Navigation Satellites) 
give the latitude dependence of electron content for essentially constant 
time (the typical duration of a high elevation NNSS pass is 20 minutes). 
Satellites in high orbits (e.g., the satellites of the Global Positioning System - 
GPS - in 20,000 km height) mix temporal and spatial dependence of 
electron content. 

The Differential Doppler data used in this study were gained by means of the 
signals of the NNSS. The receiving station was Undau/Harz, Federal 
Republic of Germany (51.62° north, 10.09° east). For a part of the satellite 
passes, the initial phase value was determined by the two stations method 
(Leitinger et al„ 1975) with data from the receiving station Graz, Austria 
(47.08° north, 15.45° east). For the remainder the initial phase value was 
determined using assumptions about the latitudinal structure of the 
ionosphere (single station method, see Leitinger and Putz, 1978). The 
latitude dependence of electron content was gained with a fixed mean 
ionospheric height of 400 km, using the geographic coordinates of the point 
in this height on the straight line from the receiver to the satellite. 

Electron content was gained from the IRI by means of numerical integration 
from 60 km to 1000 km using a step-size of 10 km. The comparison is 
based on monthly medians of electron content for the latitudes 60° north, 
55° north, 50° north, and 45° north. Because of the irregular distribution of 
NNSS passes, the day is divided into nine three-hour intervals starting at 
midnight. The monthly medians were calculated in each interval both for 
low solar activity (LSA) and for high solar activity (HSA). The criterion for 
low solar activity was a monthly mean solar sunspot number R not greater 
than 40 (nominal value R = 20), and for high solar activity an R between 130 
and 170 (nominal value 150). For LSA all months from 1975 to 1976 were 
included; for HSA a selection of months was taken from the years 1978 
through 1982 (Feichter et al., 1988, 1989). 


2.2.2 Comparison of Monthly Medians 

The comparison of electron content measurements with IRI data is 
demonstrated for a latitude of 50° north, 15° east, a typical midlatitude 
location both geographically and geomagnetically. The results of the direct 
comparison of monthly medians is shown in Figures 1 to 3. Examples for 
the seasonal variation are displayed in Figure 1 for the time intervals 06-09 
LT and 1 1-13 LT, both for low and high solar activity. Because of the quick 
change of electron content in the morning time interval, the comparison for 
06-09 LT is not too conclusive, but it can be clearly seen that the IRI gives 
an October maximum both for LSA and for HSA not present in the NNSS 
results for LSA. For noontime there are considerable differences in winter 
(Figure lc): The IRI data are too high when compared with the NNSS data. 
During summer good agreement is found for HSA (Figure Id). 
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Fig. 1 . Comparison of the annual variation of electron content (TEC) from IRI with that from NNSS 
Differential Doppler measurements in units of 10 15 mr 2 : 


a) (top left) Low solar activity, morning 
c) (bottom left) Low solar activity, noon 


b) (top right) High solar activity, morning 
d) (bottom right) High solar activity, noon 
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Fig. 2. Comparison of the diurnal variation of electron content (TEC) from IRI with that from NNSS 
Differential Doppler measurements in units of 10 15 rrr 2 : 


a) (top left) Low solar activity, January 
c) (bottom left) Low solar activity, March 


b) (top right) High solar activity, January 
d) (bottom right) High solar activity, March 






Comparing the shape of the annual variations around noon reveals that for 
HSA the absolute IRI maximum occurs in November whereas the NNSS 
maximum is seen in October (Figure Id). Since the IRI F2-peak exhibits the 
same November maximum, it Is clear that this is a feature of the CCIR 
coefficients. 

The diurnal curves exhibit the strongest differences in winter for LSA 
(Figure 2a). The shape of the winter curves is, however, similar. The 
diurnal curves for summer have different shapes (Figures 3a and b): IRI 
exhibits a double-peaked diurnal curve from May through August for LSA and 
from June through July for HSA. For low solar activity the evening maximum 
which occurs around 20 LT is even higher than the morning maximum 
occurring around 10 LT (Figure 3a). 


2.2.3 Electron Content Differences 

Since the initial value could add a bias to the NNSS data, it is advisable to 
compare electron content differences for suitably chosen latitudes. We have 
chosen the difference between the electron contents at 45° north and 60° 
north, which is comparatively free of latitudinal biases. The statistical 
properties of daytime (10-17 LT) electron content differences from NNSS 
observations are shown in Figure 4 in the form of "boxplots" (Kleiner and 
Graedel, 1980) with the range of IRI data added to the display. For LSA the 
IRI data are clearly too high, but for HSA they are too small. From October 
through February the TEC differences from IRI are even predominantly 
negative, a finding that contradicts the widespread experience that 
midlatitude TEC Increases with decreasing latitude. 


2.2.4 Equivalent Slab Thickness T and Shape Factor j \ 

Since electron content can be separated into peak electron density times 
equivalent slab thickness T and since IRI takes peak electron densities from 
the CCIR (1967) model coefficients, it makes sense to investigate the 
equivalent slab thickness T as predicted by IRI. Figure 5 shows the annual 
variation of the IRI equivalent slab thickness for different latitudes at noon 
and midnight and during low and high solar activity. For high solar activity a 
comparison was made with experimental data gained by the combination of 
electron content from the Faraday effect on the VHF signal of the 
geostationary satellite SIRIO (receiving station Graz, Austria, with the 
ionospheric point at 400 km having the coordinates 42.0° north, 10.2° east) 
with peak electron density from the ionograms of Rome, Italy (41.8° north, 
12.5° east) (Skedel, 1989). The comparison was done with monthly 
medians for the years 1981 and 1982. The results can be summarized as 
follows: 
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Fig. 4. Monthly boxplots for the difference of electron content at 45’ north and 60’ north from NNSS 
observations during daytime (10-17 LT) in units of 10 lS rrr 2 . Left hand panel: low solar activity (LSA, R = 
20). Right hand panel: high solar activity (HSA, R = 150). Boxes: from quartile to quartile, divided by the 
mecian; vertical lines to the maximum and to the minimum, respectively. Short horizontal dashes: 
uppermost and lowest deciles and confidence limits (5% error probability) of the median. The number of 
data per month is displayed at the bottom. At the left of each box, one finds the range of IRI results for the 
local time interval 1 0-1 7 LT (x — x). 


• 00 LT: The annual curve of IRI-t is below the "observed" curve. For 

1981 a decrease of IRI-t by 50 km would improve the situation 
considerably. 

• 06 LT: The IRI annual curve is still too high, but the experimental slab 
thickness is not very reliable in winter when sunrise is near 06 LT. 

• 12 LT: In summer IRI-t fits the observed t comparatively well, but in 
winter the IRI values are too high by up to 70 km. 

• 18 LT: The IRI data for 1982 were larger than the observed values 
throughout. For 1981 this is true for all months except April and June, 
which showed comparatively good agreement. The differences are 
highest in winter. 

The IRI slab thickness shows a peculiar behavior near the geomagnetic 
equator (Figure 5a) and near the crest of the equatorial anomaly (Figure 5b). 
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Fig. 5. Annual variation of equivalent slab thickness from IRI in km for different geographic latitudes and 
geographic longitude 15* east: 

£0 (top left) Vicinity of the geomagnetic equator (1 1 * north) 

b) (top right) Vicinity of the crest of the equatorial anomaly (20* north) 

c) (bottom left) Latitude 42* north (comparison with data Faraday Graz— ionosonde Rome) 

d) (bottom right) Latitude 50’ north 
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The equator T Is higher for LSA than for HSA and the equator slab thickness 
is higher than the crest slab thickness (Figures 5a and b). The crest T and 
the equator t for HSA exhibit almost no seasonal dependence. There Is 
almost no solar activity dependence of slab thickness near the crest of the 
equatorial anomaly (Figure 5b). A decrease of solar activity dependence of 
IRI-T with decreasing latitude is seen in midlatitudes, too (Figures 5c and d). 

Since the contributions to slab thickness from the electrons above the F-layer 
peak are much higher than the contributions from the sub-peak electrons, it 
is clear that the IRI deficiencies in the slab thickness behavior stem from the 
Bent model (Llewellyn and Bent, 1973). This impression is enhanced when 
one studies the behavior of other parameters of the electron density profile. 
One such parameter Is the so-called shape factor, a quantity needed for the 
assessment of higher order range errors, which play an important role in 
space geodesy and similar applications of the transionospheric propagation of 
radio waves (Leitinger and Putz, 1988). The shape factor 77 is defined by 



IRI gives a shape factor that is nearly independent of latitude and time and 
Is about 0.616 for LSA (R = 20) and 0.638 for HSA (R = 145). Since 77 
reflects strongly the electron density profile above the peak, this behavior 
means that the Bent model on which the IRI is based in this region has a 
shape that is practically independent of latitude and time and changes only 
slightly with solar activity. Recently, Buonsanto (1989) has compared 
electron density profiles measured by the incoherent scatter radar at 
Millstone Hill, Massachusetts, with IRI profiles. The comparison shows that 
even under geomagnetically quiet conditions the measured profiles can 
differ considerably from the IRI profiles, especially above the peak. Table 2 
gives the electron content TEC, slab thickness T, and shape factor 77 

obtained by numerical integration of the values read from Buonsanto' s 
figures. 


Table 2: Integral Parameters from Profiles Measured at 
Millstone Hill (46.2* N, 71.5* W) 



Low Solar Activity 

High Solar Activity 

Season LT 
HH 

Date TEC T ^ 

YYMMDD 10 15 m -2 km 

Date TEC T _ 

YYMMDD 10 15 m - 2 km 

Summer 00 
Summer 1 2 
Winter 00 
Winter 1 2 

850813 16.8 172 0.576 

850813 72.8 269 0.637 

860209 13.6 225 0.630 

860209 99.9 192 0.561 

800608 156.9 258 0.631 

800608 283.5 328 0.586 

801220 67.3 277 0.625 

801220 432.7 211 0.581 
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2.3 Morphology of Electron Temperature Anisotropy in the 
F-Region 

Koh-ichiro Oyama 

Institute of Space and Astronauttcal Science 

3-1-1, Yoshinodai, Sagamihara, Kanagawa, 229, Japan 

Abstract— Electron temperature was measured by means of planar probes 
that were Installed on Japan's ninth polar-orbiting satellite OHZORA. The 
measurements were carried out at the same time both along (7» and vertical 
(TJ to the geomagnetic field and showed the difference between these two 
temperatures. T f is more often higher than T ± during nighttime. In the 
daytime T# is sometimes lower than T± and sometimes higher. Dependence 
of the anisotropy upon geomagnetic latitude, local time, season, and height 
is investigated. 


2.3.1 Instrument and Data 

Japan's ninth scientific satellite OHZORA was launched on February 14, 

1984, to study, besides other mission goals, plasma phenomena in polar and 
equatorial regions. The polar orbiting satellite measured In the altitude 
range from 350 km to 900 km. Electron temperature probes were installed 
together with other plasma instruments, such as particle analyzers, a wave 
receiver, an electron density probe, and a topside sounder. One of the 
objectives of the electron temperature measurements was the study of the 
anisotropy of the electron temperature. Results of this investigation are 
presented here. 

Four planar probes (diameter: 120 mm) were put at the very end of the four 
solar cell paddles, as shown in Figure 1. Three probes were dedicated to 
electron temperature measurements, and one was used as a fixed bias 
Langmuir probe to measure the relative electron density. During most of the 
time, the satellite was non-spun, but occasionally it was spun very slowly. 

The maximum spin rate in the past has been 0.6 rpm. Usually, the spin rate 
was 0.2 rpm or less. The satellite spin axis was controlled to point towards 
the sun with an accuracy better than 1°. 

The principle of the instrument and the data reduction method was 
described by Hirao and Oyama (1970) and Oyama et al. (1985). The probes 
whose sensor plates are vertical to the geomagnetic field lines (S3 and S7 in 
Figure 1) measure the parallel electron temperature Tf , and the probes 
whose sensor plates are parallel to the geomagnetic field lines (SI and S5 in 
Figure 1) measure the vertical electron temperature T±. For this statistical 
study only those data points were used for which at least one of the sensors 
was aligned with the magnetic field to within ± 20°. The data were also 
carefully checked to remove measurements carried out in disturbed areas in 
the satellite wake or in the thick sheath which may be caused by the strong 
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Fig. 1. Satellite OHZORA and location of the four temperature probes; sensors Si, S3, and S5 are used 
as the electron temperature probes while the sensor S7 measures ion and electron currents. 
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satellite charging. Altogether, 505 data points are used from observations 
carried out between February 1984 and December 1985. 


2.3.2 Results for Individual Orbits 


Figure 2 shows results which were obtained with sensors SI and S3 during 
the satellite orbit 947; the two sensors are orthogonal to each other. The 
upper panel shows the angles between the geomagnetic field direction and 
the three satellite-centered axes X, Y, and Z (see Figure 1). To simplify 
matters, we will call these angles X, Y, and Z as indicated in Figure 2. Also 
shown (thin line) is the floating potential, which is measured by means of 
one of the sensor electrodes with respect to the satellite potential. The 
second panel from the top shows the ratio of electron temperatures 
obtained by the two sensors SI and S3. The individual electron 
temperatures are plotted in the two lowest panels. Orbit attitude 
information is listed at the bottom of Figure 2, Including the Universal Time 
(UT), the altitude (HHGT), L-value (FL), geomagnetic latitude (GLAT), 
geographic latitude (XXLAT), longitude (XXLON), and magnetic field 
strength (BB). 

From 0051 UT to 0056 UT, the satellite traversed the polar cap, observing a 
large scatter of electron temperature values during that time. Oyama and 
Abe (1987) concluded that this scatter is real and reflects the geophysical 
conditions in the polar ionosphere. Between 0056 UT and 0057 UT, both 
temperatures exhibit a dip decreasing to a minimum value of 2000 K at 
005640 UT and Increasing again to 2600 K at 0057 UT. This region 
corresponds to the boundary where the flux of precipitating auroral 
electrons (0.2~16 keV) starts decreasing and the flux of trapped energetic 
particles of 0.1-140 MeV energy begins to increase. Simultaneous 
measurements show a maximum of the counting rate of energetic particles 
at 0058 UT. & ^ 

At the beginning of the orbit the angle Yis about 170°. This means that 
sensor S3 is nearly vertical to the geomagnetic field, and, therefore, it is 
picking up most of the parallel component of electron temperature T§, 
while the sensor SI is picking up the vertical component T±. The panel 
second from the top shows that T(S1)/T(S3) is 0.7 at this time; that is, T, is 
1.4 times higher than T±. From 0051 UT to 0059 UT, the ratio increases 
toward one as the angle changes from 170° to smaller values. At 0059 UT 
the angle Y is nearly 135° and the angle Z is 108°. In this configuration SI 
and S3 should pick up equal fractions of the electron temperatures T# and 
Tj_. In fact, the second panel shows that T(S1)/T(S3) is nearly equal to one 
at this time. After 0059 UT the situation reverses, and the parallel 
temperature is more strongly observed by sensor SI. T[S1)/T(S3) is 1.1 at 
0102 UT indicating that the parallel component is still larger than the 
vertical component. 
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Fig. 2. Electron temperatures measured by the planar probes SI and S3 during the pass 947 (lower two 
panels). Top panel shows the angles X, Y, and Zbetween the X, Y, and Z axes (see Figure 1 , page 24) 
and the geomagnetic field direction. Also shown is the floating potential. The second panel shows the 
ratio between the two temperatures measured by sensors Si and S3. Numbers at the bottom are 
Universal Time in hours and minutes (UTIME), height in km (HHGT), L value (FL), geomagnetic latitude 
(GLAT), geographic latitude and longitude (XXLAT and XXLON), intensity of geomagnetic field in Gauss 
(BB), and IW (sunlight indicator: IW = 0 means satellite is sunlit). 
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Simultaneous measurements of electron flux In the energy range 0.2 to 16 
keV show a strong precipitation pattern during the polar cap crossing from 
0051 UT to 0057 UT. Large ion counting rates are seen from 0051 UT to 
0056 UT. Thus, there is a time period and region in which electrons are 
detected, but no ions. This period from 0056 UT to 0057 UT corresponds 
to a dip in electron temperatures and a slightly enhanced electron 
temperature anisotropy. 

Figure 3 shows the temperatures and the anisotropy observed during the 
satellite pass 2853. The angles Zand X are about 100° and 15°, 
respectively, at the beginning of the pass. In this configuration, sensor SI 
detects most of the Tf component, while S3 detects most of the T± 
component. The ratio of these two temperatures Is about 1.5; that is, T# is 
about 1.5 times higher than T ± . Between 2336 UT and 2337 UT, the 
configuration is such that both sensors should see about equal fractions of 
the two temperature components. Figure 3 shows that, indeed, the ratio of 
the electron temperatures from the two sensors is nearly one during this 
time period. 

For our statistical analysis of the electron temperature anisotropy, we use 
only electron temperatures obtained when the normal of one planar 
electrode deviates less than 20° from the magnetic field direction and when 
the vertical component can be simultaneously obtained from an orthogonally 
mounted sensor. Figure 4 shows T# versus T± for the orbit period from 
2331 UT to 2333 UT in Figure 3. During this time period, the angle 
between the normal of sensor SI and the geomagnetic field changes from 
18° to 24°, and, therefore, sensor SI detects mostly T§ while sensor S3 
detects mostly T±. The data can be approximated by the relation T| = 1.2 
T± as shown by the solid line in Figure 4. Figure 5 shows data taken 
between 0116 UT and 0118 UT during the satellite orbit 2721. The broken 
line shows the best fit curve, which is T, = T ± + 250. In this case, the 
average temperature ratio (T f = 1.07 T ± solid line) is obtained by giving 
more weight to the data points for which the sensor SI is aligned the 
closest to the geomagnetic field. 


2.3.3 Statistical Results for the Occurrence of Anisotropy 

From figures similar to Figures 5 and 6, we have calculated the temperature 
ratio R (= Tf /T ± ) for 505 satellite passes and investigated the dependence 
of R upon geomagnetic latitude, local time, season, and altitude. Figure 6 
shows histograms of R for different geomagnetic latitude ranges. All data are 
sampled into six geomagnetic latitude zones without discriminating 
Southern and Northern Hemispheres. The average height for each data 
group varies from 555 km to 470 km. At latitudes less than 20°, the 
number of cases with R less than one is almost equal to the number of cases 
^th R greater than one. The number of cases with R greater than one 
increases as we go to higher latitudes. We conclude that anisotropy is 
becoming more likely as we go to higher latitudes. 
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Fig. 3. Same as Figure 2 for orbit 2853. 
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REV. 2853 8/ 23 233100- 233300 



Fig. 4. Parallel electron temperature versus vertical electron temperature during orbit 2853 and best 
linear fit (7# = 1.21 TJ. 


REV. 2721 8/ 15 U6Q0 - 11800 



Fig. 5. Parallel electron temperature versus vertical electron temperature during orbit 2721 and best 
linear fit 7> =T ± + 250 (dotted line); also shown is 7>= 1.07 Tj_ (thick line). 
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Secondly, all data were distributed into six local time ranges. Each of the 
six data sets consists of measurements taken during a four-hour period as 
shown In Figure 7. The average height ranged from 515 km to 646 km. 
Each of the six cases indicates a temperature anisotropy with R greater than 
one, that Is, T # greater than T±. The temperature ratio seems to Increase 
towards midnight. 

In order to study how the anisotropy varies with height, all data are grouped 
into six height zones starting from 300 km In steps of 100 km, as shown in 
Figure 8. Although we cannot see the influence of height upon the electron 
temperature anisotropy clearly, a weak tendency seems to exist that T# 
approaches T± gradually as the height increases. It should, however, be 
noted that temperature measurements become more difficult with 
increasing height because of decreasing electron density. 

Finally, Figure 9 shows the annual variation of R. All data are grouped into 
six bimonthly periods, starting January 1985. These histograms seem to 
show that the anisotropy T# > T± appears more often in summer (July and 
August) than in winter (January and February). 

Similar statistical analysis is planned to examine the dependence of JR on 
longitude and on solar and magnetic activity. 


2.3.4 Conclusions 

The investigation of electron temperature anisotropy observed by the 
OHZORA satellite can be summarized as follows. 

• In the equatorial region, Tj is sometimes higher than T± and sometimes 
lower. With Increasing latitude, the anisotropy T # > T ± dominates. 

• T f Is higher than T± for all local times. The difference increases and 
becomes clearer toward evening and during nighttime. 

• Anisotropy (Tf /T ± >1) appears more often in summer than in winter. 

These results clearly show that anisotropies have to be considered In 
describing the global and temporal variations of the ionospheric electron 
temperature. It is therefore recommended that future satellite instruments 
use temperature probes that can discriminate between the two temperature 
components. 
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Fig. 6. Histograms of R = 7j / 7j.for different latitude zones. 
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Fig. 7. Histograms of R = 7ji / T± for different local time periods. 
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Fig. 8. Histograms of R = Tg / T±\or different height ranges. 
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2.4 Model Descriptions for the Ion Transition Heights 
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I. Kutiev 
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Abstract— The ion composition in the upper and middle ionosphere exhibits 
two characteristic transition heights that mark the upper and lower boundaries 
of the region dominated by atomic oxygen ions. The upper transition height Hj 
is found at the point where the ion gas consists of equal parts of oxygen and 
light ions, and the lower transition height ht where the oxygen density equals 
that of the molecular ions. We investigate the global and temporal morphology 
of these heights, summarizing the results of earlier studies. The use of these 
characteristic points for the IRI ion composition model is emphasized. 


2.4.1 The Upper Transition Height 

The ion gas of the topside ionosphere consists mainly of 0 + , H + , and He + ions. 
0 + and He + ions are produced by photoionization, whereas the main source of 
H + ions is the charge exchange reaction between H and 0 + . Production and 
loss together with redistribution by diffusion and ion drift are the major 
processes that shape the distribution of the different ion species. Figure 1, 
taken from the work by Raitt et al. (1975), shows 0 + and H + profiles for 
different upper boundaiy conditions. As a result of the difference in ion mass 
and because of the charge exchange reaction, the amount of 0 + ions decreases 
with altitude, whereas the percentage of light ions increases. The transition 
height Ht is defined as the height at which equal parts of 0 + and light ions are 
present. Ht depends strongly on the ion flow conditions at the plasmaspheric 
boundaiy, as is shown in Figure 1. Thus, Ht is a critical parameter in specifying 
and assessing the topside ion gas. Several studies (see Table 1) have examined 
the global and temporal variation of this characteristic height. 

In summary, these investigations established the following characteristic 
variation patterns: 

• The diurnal variation of Ht follows the change in 0 + density, being lowest 
at night and largest during daytime. Ht can drop down to as low as 500 
km during nighttime and is rarely below 800 km during daytime. 

• During daytime Ht is lowest at the magnetic equator and increases 
towards higher latitudes. This is caused by a similar change in plasma 
temperatures. Higher temperatures shift the ion density profiles upward. 
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ALTITUDE (KMT 



Fig. 1. Theoretical H+ and 0 + density profiles as a function of altitude 
km. The H + velocities at the upper boundary in units of 

3 0 m 5 0 (at 10 0 and (h) 20.0. The shaded region shows the range of 0 + densities tne lower 
^unda^wrresporK^ng to H- profile (a) and the upper boundary to (h) (Raitt et al. ( 1975). 
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Table 1: H T Studies and Typical Values 


Author 

Data 

Used 

■I 

[ ~ H T 7km — i 

Day 

Night 

Geomagnetic Latitude 
Low Middle 

Geomagnetic Latitude 
Low Middle 

Brinton 
et al. (1970) 

Explorer 32 
IMSt 

Low 

1300 

1900 

- 

- 

Goel et al. 
(1970) 

OGO 4 
Explorer 31 
RPA 

Medium 

1300 

1600 

700 

700 

Titheridge 

(1976) 

Alouette 1 
ISIS 1 
Topside 
Sounder 

Low 

800 

(800) 

1200 

(1200) 

650 

700 

(550) 


m« 

Low 



■kI 

550 

(550) 

Kutiev et al. 
(1984) 

OGO 6 
RPA 

mm 



1100 

800 

Bilitza (1984) 

Incoherent 

Scatter 

Low 

1400 

>1500 

700 

1100 


Summer value; winter value in brackets. 
■’ IMS - ion mass spectrometer, 
tt RPA - retarding potential analyzer. 


• Ht is higher in summer than in winter by 100 km to 200 km. An 
interhemispherical proton flux from the summer to the winter 
ionosphere is thought to be the reason for this difference (Brinton et al., 
1970). 

• Hj increases with solar activity. Titheridge (1976) finds nighttime 
transition heights about 500 km during solar minimum and about 1100 
km at solar maximum. 

There are, however, still quite a number of open questions and conflicting 
results that have to be resolved by future studies and measurements, for 
example: 

• Table 1 shows that the earlier satellite measurements indicated 
generally higher transition heights than the more recent ones. There is 
also an obvious discrepancy between transition heights obtained by the 
incoherent scatter technique and those obtained by the in situ 
techniques. 

• Kutiev et al. (1984) report considerable longitudinal changes in 
nighttime Ht obtained from OGO 6 RPA measurements, whereas the 
Taiyo RPA data (Miyazaki, 1979) exhibit a negligible longitudinal scatter. 


35 


































Only a few investigators have attempted to represent their results by simple 
empirical relationships. Miyazaki (1979) used the following formulas to 
describe transition heights deduced from the Taiyo RPA measurements 
between March 1975 and October 1976. 


Hr/km = j 


591 -87.1 sin 2 T 

.591 - 87.1 sin 2 T + (217 + 473 sin 2 ¥)(cos X) 0 - 331 


X > 90° 
else 


T is the magnetic dip latitude and X the solar zenith angle. Kutiev et al. 
(1984b) used a power series in longitude and in magnetic dip latitude to 
represent the nighttime H T deduced from OGO 6 RPA data. 

In Figure 2 we compare the transition height predicted by these two H T 
models with transition heights obtained from the IRI ion composition 
models and from the ion composition model by Danilov and Yaichnikov 
(1985). For more information about the latter two models, see section 3.4, 
page 71, of this document. The diurnal variation of Hr is shown at the 
magnetic equator (Lima) and at a magnetic midlatitude location (Arecibo). 
The IRI model predicts considerably lower transition heights than the other 
models. It also produces unreasonably abrupt changes of H T at sunrise and 
sunset. The nighttime H T from the model by Kutiev et al. (1984b) is much 
higher than all the other predictions. Ht deduced from the ion composition 
model by Danilov and Yaichnikov (1985) agrees fairly well with Miyazaki's 
(1979) H T model at midlatitudes. At the magnetic equator the need for 
further data acquisition and modeling work is most obvious: None of the 
four models has much in common with the others. 


2.4.2 The Lower Transition Height 

In the altitude range 100 km to 300 km the ionospheric ion gas consists 
mainly of 0 + , 02 + , and NO + ions. In the absence of electric fields the 
distribution of ions can be explained straightforwardly by the chemical 
balance of photoionization and recombination. As a result of these 
processes, the concentration of 0 + ions increases with altitude whereas the 
percentage of molecular ions decreases. At the transition height ht the 
percentage of molecular ions has dropped to 50%. 

Oliver (1975) investigated the diurnal variation of this characteristic height 
by examining 40 rocket IMS measurements from 1956 to 1971. He found 
that the data could be well represented by 

h t / km = hs + 203.7 + 24.7 tanh(- (88.4 - *)/19.7) 

The remaining seasonal variation was approximated in terms of annual and 
semi-annual terms varying with the day number d. 

hs / km = - 16.8 cos((d + 22) 2n/365) - 9.4 sin((2d + 10.5) 2 tc/365) 
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local time 


Fig. 2 a, b. Upper transition heights H T obtained from the models by Miyazaki (1979), Danilov and 
Yaichnikov (1985), Kutiev et al. (1984), and from the IRI at the magnetic equator (top) and at a magnetic 
midlatitude location (bottom). Also shown is the solar zenith angle (plus 300). 
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Lathuillere and Brekke (1985) report good agreement between their 
EISCAT incoherent scatter results and Oliver's formula. 

In Figure 3 the ht values obtained with Oliver's (1975) model are compared 
with transition heights predicted by the IRI and by the ion composition 
model of Danilov and Yaichnikov (1985). At midlatitudes IRI agrees fairly 
well with Oliver's results, whereas the model by Danilov and Yaichnikov 
(1985) exhibits an unreasonably small change from day to night. At the 
equator IRI predicts the transition heights ht much lower than the other 
two models. 


2.4.3 Conclusion 

The lower and upper ion transition heights are parameters which 
characterize and specify the ion composition profiles throughout most of the 
ionosphere. Our review has found only a few studies which examined the 
global morphology and the temporal change of these heights. Unfortunately, 
these studies have revealed somewhat conflicting results, in particular close 
to the magnetic equator. All three techniques used to evaluate transition 
heights have their specific accuracy problems. The Ion mass spectrometers 
(IMS) and retarding potential analyzers (RPA) flown on long-lasting satellites 
missions have to deal with calibration problems and in-flight contamination 
effects. Ground-based incoherent scatter radars need long integration times 
or some assumptions about the plasma temperatures to deduce ion density 
ratios from the backscattered signal (Lathuillere et al., 1983). 

Inconsistencies can only be resolved by comprehensive statistical studies of 
the data base accumulated so far and by a continued measurement program 
for these parameters from ground and space. It is clear that a better 
understanding of the variation patterns of the transition heights will 
significantly improve the modeling and prediction of ionospheric ion 
composition. 


38 



I.IMA r/vw i:>s1 siliuukit i?-^ i r»r» 



local lime 



local time 


Fig. 3 a, b. Lower transition height fy obtained from the models by Oliver (1975), Danilov and 
Yaichnikov (1985), and from the IRI at the magnetic equator (top) and at a midlatitude location (bottom) 
Also shown is the solar zenith angle (divided by 1 0; plus 150). 
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2.5 Models for Horizontal E- and F-Region Drifts 


E. S. Kaztmirovsky, E. I. Zhovty, and M. A. Chernigovskaya 
Siberian Institute of Terrestrial Magnetism, Ionosphere, 
and Radio Propagation (SibIZMIR) 

664033, Post Box 4026 
Irkutsk, U.S.S.R, 

Abstract — E- and F-region horizontal ionospheric drift data, obtained 
worldwide from 1957 to 1970 by ground-based radio methods D1 and D3, 
are statistically analyzed and a computer code for the average variations in 
latitude and local time is constructed. This model has been adopted for the 
International Reference Ionosphere, IRI-90. 


2.5.1 Data and Method 

Plasma transport is very important for understanding the variations of the 
ionosphere in space and time. Therefore, the URSI/COSPAR Task Group on 
IRI has initiated an effort to construct a computer code that generates ion 
drifts for chosen location and local time. 

A large quantity of experimental data obtained by measurement of horizontal 
ionospheric drifts is now available on a worldwide scale. The data we 
considered here were obtained by the spaced -receiver methods D1 and D3, 
as described by Briggs (1977). We have developed appropriate procedures 
that allow us to infer from these data the main parameters of the global 
ionospheric motion at E- and F-region levels (Vergasova et al., 1979). 

Thus, following our earlier proposals (Kazimirovsky et al., 1984, 1985), we 
present here an analytical description of the dynamic behavior of the 
terrestrial ionosphere. 

The experimental data are from 23 stations in the Northern Hemisphere 
and were obtained between 1957 and 1970. The worldwide coverage in 
geographic latitude is 7° to 71° north (7.5° to 64.1° north in geomagnetic 
latitude) and 0° to 131° east in geographic longitude. 


2.5.2 The Drift Models 

About 100000 individual drift velocity values in the E-region (90 km < h < 
135 km) and F-region (h > 135 km) were statistically treated with sampling 
according to season and to solar activity as indicated in Table 1 and Table 2. 
Definition of seasons is as follows: December through February = winter, 
March through May = spring, June through August = summer, September 
through November = fall. 
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"SEASON = all" corresponds to the averaged velocities for the whole 
observed period. Three periods of solar activity are considered: (1) the 
whole solar activity range from minimum to maximum (all), (2) solar 
maximum [FI 0.7 = 200), and (3) solar minimum [FI 0.7 =70). FI 0.7 is the 
solar radio flux at 10.7 cm wavelength (Covington index). 


Table 1: Analysis Periods for E-Reglon 


PERIOD 

1 

2 

3 

4 

5 

6 


winter 

sprinq 

summer 

fall 

winter 

spring 

F10.7 

all 

all 

an 

aU 

200 

200 


PERIOD 

7 

8 

9 

10 

11 

12 

SEASON 

summer 

fall 

winter 

sprinq 

summer 

fall 

F10.7 

200 

200 

70 

70 

70 

70 




Table 2: 

Analysis Periods for F-Reglon 



PERIOD 

1 

2 

3 

4 

5 

6 

7 

SEASON 

- 

winter 

sDrina summer 

fall 

winter 

sprinq 

F10.7 

all 

all 

all 

all 

al 

200 

200 


PERIOD 

8 

9 

10 


11 

12 

13 

SEASON 

summer 

fall 

winter 


sprinq 

summer 

fall 

F10.7 

200 




70 

70 



An algorithm for computing the zonal and meridional drift components VX, 
VY can be found in section 3.5, page 75, of this document. It is written in 
PASCAL 8000 for IBM 360/370 compatible computers. The main block of 
the program is the procedure DRIRR for calculating VX and VY for a certain 
period (P), geomagnetic latitude (FI), and local time (LTT). VX and VY are 
first computed by Fourier development: 

3 

VX (LTT, FI) = Xo(FI) +£ X n (FI) sin(n LTT + CX n (FI)) 

n=l 

3 

VY (LTT, FI) = Y 0 (FI) + Y Y n (FI) sin(n LTT + CY n (FI)) 


for six main latitudes FI = 10°, 20° 60°. Spline interpolation between 

these allows determination for any specified value of FI. The coefficients 
X n (FI), Y n (FI), CXn(FI), and CY n (FI) are contained in the files FAR1 and FAR2 
listed in section 3.5, page 75, of this document. 
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3.1 Introduction 


IRI-90 describes monthly averages of electron density, electron 
temperature, ion temperature, and ion composition in the non-auroral, quiet 
ionosphere from 100 km to 1000 km. 


3.1.1 Data Sources 

The following data sources were used to establish IRI-90. 


• Worldwide lonosonde measurements (peak plasma frequencies => peak 
densities) 

• Incoherent scatter measurements from Jicamarca, Arecibo, St. Santin, 
Millstone Hill, and Malvern (E-valley, F2-peak height, F2-bottomside, 
topside, electron and ion temperature) 

• Alouette 1, 2 topside sounder measurements (topside) 

• AE-C, ISIS 1 and 2 In situ measurements (electron temperature) 

• AEROS in situ measurements (electron and ion temperature. Ion 
composition) 

• Rocket measurements (D- and E-region electron density and ion 
composition) 

• Beacon satellite data (ionospheric electron content for test purposes) 

® Ground-based absorption measurements (D- and E-region variability) 


3.1.2 Functions 

A variety of analytical expressions and functions is used in IRI-90 to 
represent the temporal and spatial variations of the ionospheric densities 
and temperatures. Global variations (with latitude and longitude/local time) 
are in most cases described by a form of a spherical harmonics series 
(Legendre polynomials). 

A j?J° jP °f f unctions introduced by S. Epstein has been used to represent 
altitudinal variations as well as special latitudinal and diurnal features. The 

nrst three members of the Epstein family of functions are shown in Figure 1 
and defined as & 


EPS_! (h ; HX , SC ) = ln(l + e*) 
EPS 0 (h ; HX, SC ) = 1/(1 + e -*) 
EPSi (h ; HX. SC ) = eV(l + e *)* 


(3.1a) 

(3.1b) 

(3.1c) 
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Fig. 1. Epstein functions. 


with 


x = (h - HX)/ SC 

(3. Id) 

EPSi+i = d EPSi/dx 

(3. le) 


EPS_i describes a transition, EPS 0 a step, and EPSi a peak at h = HX with 
the width SC. Three combinations of Epstein functions are explained in the 
following sections. All three have been used repeatedly in constructing the 
IRI profiles. 

The Epstein step function can also provide a type of filter function taking 
the value 0 below a certain altitude and the value 1 above. Rawer (1987) 
explained how three such filter functions can be used for a joint analytical 
representation of the whole electron density profile. 

3. 1.2.1 Booker Profile Function 

A parameter Y, which steps from a value Yi to a value Y 2 at h = hs, can be 
described by the function 
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EPSTEP(h; Y lt Y 2 , hs, d s ) = Yj + (Y 2 - Yj) EPSo(h; hs, d s ) 


(3.2) 


d s determines the width of the transition region. Booker's (1977) approach 
uses several such steps. He assumes that the altitudinal profile of a 
parameter y can be divided into several subsections such that in each 
subsection y can be fairly well approximated by a straight line or, in other 
words, that the gradient dy/dh is nearly constant in every subsection (e.g., 
Figure 5). Booker calls this the skeleton profile. The derivative of the 
skeleton function can then be represented by a sum of EPSTEP functions, 
stepping from subsection to subsection. Integrating from ho to h, we obtain 
y in terms of the "Booker function" B. 

y(h) - y 0 = B(h; ho. M, DY, HX, SC) 

M 

= (h - ho) DYi + (DY J+1 - DYj) SCj ( EPS^fh; HXj, SCj) 

- EPS_i (ho; HXj, SCj)) (3.3) 

M is the number of subsections, DYj [j = 1 to Af + 1) are the constant 
gradients in these subsections, HXj [j = 1 to Af) are the subsection 
boundaries, and SCj (/ = 1 to M) are the thicknesses of the transition zones 
between subsections, yo is the integration constant, with y (ho) = yo- 
Suitable transition thicknesses SCj have to be found by trial and error. Small 
SCj values produce a profile closer to the skeleton than to the real profile; 
large SCj values provide greater smoothness but could cause interference 
between adjacent step functions. Best results are obtained with 

D/20 £ SCj ZD/ 10 (3.4a) 

where D is the height range of the smaller of the two adjacent subsections 

D = min ( hj - hj- lt h J+1 - hj) (3.4b) 


Care has also to be taken in choosing the anchor point (integration 
boundary) ho- It should not be located in one of the highly variable transition 
regions. The best choice is a ho from a region where the skeleton and the 
original profile are close together. These constraints make it difficult to 
represent layer profiles with Booker's method. A profile peak can be 
approximate by either two skeleton lines whose gradients have opposite sign 
or by three skeleton lines, where the gradient of the middle one is zero. In 
the first case, the tip of the layer is in the middle of the transition zone and 
the peak height cannot be used as anchor point ho . In the second case, 
condition (3.4) will almost always be violated. 
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3. 1.2.2 Rawer Lavfer) Function 

A better way of representing layer shape was introduced by Rawer (1984). 

His LAY function combines the Epstein transition function with a linear term. 

LAY (h; hm, HX, SC) = EPS_i(h; HX, SC) - EPS-i (hm; HX, SC) 

- (h - hm) EPSo (hm; HX, SC) /SC (3.5) 

It peaks at h = hm, where it takes zero value. A wide variety of layer shapes 
can be represented by the LAY function as shown in Figure 2. The 
appropriate HXand SC values can be found by a non-linear least-square-fit 
procedure. 

3. 1 .2.3 Dav-Night Transition Function 

Epstein step functions are also used to describe simple day-night variations 
of ionospheric parameters. The combination of two step functions 

HPOL(t; Y d , Y n , SR, SS) = Y N + [Y D - Y N ) EPSqU; SR, TD) 

+ [Yn - Y d ) EPS off, SS, TD) (3.6) 

varies with local time t continuously from a constant nighttime value Y N to a 
constant daytime value Y D . The steps occur at the local times of sunrise SR 
and sunset SS and the step width is determined by TD (usually one hour). 
This local time variation is used for all parameters for which only day and 
night values are listed. 
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Fig. 2. Rawer's LAY function for a variety of HX and SC values. 





3.2 Electron Density 

The IRI electron density profile is divided into six subregions as shown in 
Figure 3, including the topside, the F2-bottomside, the FI -layer, the 
intermediate region, the E-valley, and the E-bottomside and D-region. The 
boundaries of these subsections are marked by several characteristic profile 
points including the F2-, F1-, and E-layer peaks. The critical frequencies 
foF2, foFl , and foE of these peaks have been monitored by the worldwide 
network of ionosondes since the forties. A radio wave with the critical 
frequency JoF2 transmitted vertically from the ground will be reflected at 
the F2-peak. The squares of the critical frequencies are proportional to the 
electron density at the height where reflection occurs, 

NmF2/m~ 3 = 1.24 10 10 [foF2/MHz ) 2 (3.7) 



Fig. 3. Buildup of IRI electron density profile. 


51 


preceding 


Rage blamk 


,V0T FUMED 




3.2.1 Topside (hmF2 to 1000 km) 

The topside profile is normalized to the F2-peak density NmF2 and altitude 
hmF2. 


3-2.1 .1 F2-Peak Density (AftnF2. ft?F2) 

IRI-90 offers two choices for the description of the F2 critical frequency 
foF2 (and therefore NmF2 ): (1) the model recommended by the 

International Radio Consultative Committee (CCIR, 1967) and (2) the one 
recently proposed by the International Union of Radio Science (URSI, Rush 
et al., 1989). The mathematical description is the same in both cases. It is 
based on first Fourier analysis (seventh order, 15 coefficients) of the 
monthly median diurnal variation as observed by the worldwide network of 
ionosondes, and secondly, a worldwide description in terms of spherical 
(Legendre) functions, separately found for each of the 15 coefficients (Jones 
and Gallet, 1965). Two sets of 988 coefficients each are given for each 
month, one for low and one for high solar activity (Ri 2 = 0, 100); foF2 for 
intermediate activities is found by linear interpolation. 

For very high solar activities foF2 observations indicate a saturation effect. 

In IRI this effect is enforced by keeping JoF2 constant above R 12 = 150. 

(See Equation [3.13d]). 

Because of the strong magnetic control of F-region processes, it is 
advantageous to use a magnetic field coordinate in the global analysis. The 
best results were obtained with the modified dip latitude (short: Modip), n, 
introduced by Rawer (1963) 

tan n = /cos 1 / 2 ^ (3.8) 

which at low latitudes is near the magnetic inclination (short: dip), i|), and 
gets closer to the geodetic latitude cp as latitude increases. 

For the newer URSI model, the data sparse regions (oceans. Southern 
Hemisphere) were first populated with JoF2 values obtained by aeronomic 
theory and then the analysis procedure was applied (Rush et al., 1984, 

1985). The URSI coefficient set provides a somewhat better representation 
of ocean JoF2 values than the CCIR model, when compared with ISS-b 
topside sounder measurements (Rush et al., 1989). Recent reviews of JoF2 
modeling were given by Bilitza et al. (1987) and Bradley (1990). 

3-2. 1.2 F2-Peak Height (hmF2. M(3000)F2) 

The F2-peak height hmF2 is obtained by its close correlation with the 
propagation parameter M(3000)F2 (Shimazaki, 1955; for review, see Bradley 
and Dudeney, 1973, and Bilitza et al., 1979). M(3000)F2 is defined as 

M(3000)F2 = MUF/foF2 (3.9) 
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where MUF Is the maximum usable frequency that, refracted in the 
Ionosphere, can be received at a distance of 3000 km. This factor has been 
routinely scaled from ionograms, and numerical maps (CCIR, 1967) have 
been established in the same way as described above for JoF2. hmF2 is 
calculated from M(3000)F2 with the empirical formula (Bilitza et al., 1979) 


hmF2 / km = 1 490/ [M(3000)F2 + DM) - 176 (3.10a) 

with the correction factor 

DM = fif; 2 /{foF2/foE - f 3 ) + f 4 (3. 10b) 

and the solar activity functions 

fi = 0.00232 R 12 + 0.222 (3. 10c) 

f 2 = 1 -R 22 /150 exp (-(¥ / 40) 2 ) (3.10d) 

f 3 = 1.2 - 0.0116 exp(Rj2/41.84) (3.10e) 

f 4 = 0.096 {R 12 - 25)/ 150 (3.10f) 

R 12 Is the 12-months-running mean of solar sunspot number and T is the 
magnetic dip latitude 

tanT = ^ tanijj (3.11) 


which is related to the magnetic inclination (short: dip) iJj of the Earth's 
magnetic field at 300 km altitude. 

Formula (3.10) is based on the work by Shimazaki (1955), who was the first 
to describe the strong anti-correlation between hmF2 and M(3000)F2, and 
the work by Bradley and Dudeney (1973), who were the first to take account 
of the effect of the E-layer ionization on the correlation. Bilitza et al. (1979) 
used incoherent scatter data to improve the earlier relationships, which 
were based on the, generally, less reliable hmF2 values deduced from 
ionosonde recordings. 

3.2. 1.3 Topside Profile Shape 

The topside model is based on the descriptive compilation of Alouette 
topside sounder data by Bent and his colleagues (Llewellyn and Bent, 1973). 
Epstein functions are used to obtain an analytical representation of the 
values given by Bent in tabular form. 

N (h)/NmF2 = exp(-l/a • B(x; xq, 2, TG, TX, TC)) (3.12a) 
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B is the Booker function explained in section 3. 1.2.1 (Equation [3.31) and x 
is a modified altitude variable transforming the F2-peak to near 300 km 

x - xo = a (h - hmF2 ) (3. 12b) 

with the transformation factor 

a = 700/(1000 - hmF2) (3.12c) 

and the new peak height 

x 0 = 300 - 6 (3.1 2d) 

The Epstein parameters for the topside profile are given in Table 1. 


Table 1: Epstein Parameter for Topside Formula 


/ 

TG 

TX 

TC 

i 

e 

300 

100 

2 

0 

394.5 

p 

3 

^ 

— 

— 


The parameters i, tj, p are functions of the geomagnetic latitude 4>, the 
monthly solar radio flux FI 0.7, and the F2 critical frequency JoF2 

i, q, p = t 0 + ti T m (4>) + t 2 R(F10.7) + t 3 T(4>) R(FJ0.7) + 


where 


UfoF2 + t 5 /oF2 T(4>) + t&foF2^ 

(3.13a) 

T(4>) = cos 2 4> 

(3.13b) 

fEPSi(4>, 0. 15) for q 1 
Tml4>J " Icos 2 * for i pj 

(3.13c) 

nrrrrn n _ J(FI 0.7 - 40) /30 for FI 0.7 < 193} 
R(FJ0.7)-| 51 else | 

(3.13d) 


The coefficients t t were determined by fitting these functions to the 
corresponding values obtained from Bent's tabulation and are listed in 
Table 2. 
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Table 2: Topside Profile Parameters 



11 

£ 

P 

la 

0.058798 

0.078922 

! -128.03 

ti 

-0.08 

-0.0046702 

20.253 

^2 

-0.014065 

-0.019132 

-8.0755 

^3 

0.0069724 

0.0076545 

-0.65896 

U 

0.0024287 

0.0032513 

0.44041 

15 

0.004281 

0.006029 

0.71458 


-0.0001528 

-0.00020872 

-0.042966 


The corrective term 6 Is related to rj, 4, and p by 


6 




T1 

1 + Z 


1 

2 




Z 

(1+Z)2 + 



with 


(3.14a) 


Z = exp (94.45/p) (3.14b) 

It should be noted that the parameters TG, TX, and TC for the Booker 
function were obtained by a least-square-fit procedure, rather than from a 
skeleton profile as originally intended by Booker (1977) (see section 3. 1.2.1, 
page 46). Thus, the topside representation is closer to the LAY formalism 
(see section 3. 1.2.2, page 48) as pointed out by Rawer (1987). Rearranging 
formula (3.12a) in terms of LAY functions (Equation [3.5]), one obtains 


N(h)/IVmF2 = exp (- pri/oc LAY (x; x 0 , 394.5. P) - 

100 ^ LAY (x; x 0 , 300, 100) - (h - hmF2 ) Z) (3. 15a) 


where 


Z = r|/(l + exp((6 + 94.5)/p)) + 4/(1 + exp(6/100)) - 4 (3.15b) 

In general, Z is so small that the linear term can be neglected in computing 
the argument of the exponential function, at least in the lower topside. 

Good agreement was found in comparisons of IRI with AEROS and AE-C 
satellite data and with Jicamarca incoherent scatter measurements (Bilitza, 
1985a). Buonsanto (1989) reported that IRI overestimates the electron 
densities observed by the Millstone Hill incoherent scatter radar. A 
comparison with DE 2 satellite data indicated that IRI topside density may 
be too low for very high solar activities (Bilitza et al., 1987). 
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3.2.2 Bottomside (hmFl to hmF2) 

The bottomside profile description was established by Ramakrishnan and 
Rawer (1972). It is defined as (see also Figure 4) 


N {h)/NmF2 = expt-x^J/coshU) 


(3.16a) 


with 


x = ( hmF2 - h)/B 0 


(3.16b) 


NmF2 and hmF2 Eire the electron density and height of the F2 peak, and Bj 
is equal to 3 in most cases (see section 3.2.4, page 59. for exceptions). IRI- 
90 offers two choices for the bottomside thickness parameter Bo: (1) a table 
of values (Table 3) based on reduced ionosonde profiles as in previous __ 
editions of IRI, and (2) Gulyaeva's (1987) model for the half density height 
ho 5 . In most cases, the newer model (option [2]) is the better choice, 

especially at low latitudes. 

Table 3: Bottomside Thickness Parameters B 0 Deduced from lonograms Recorded at 

Mexico City, Huancayo (Peru), and Lindau (F.R.G.) 


BqI km 


Modip 

=18 




R= 10 


f ?=100 


R= 10 


R= 100 


Spring 


LT=1 2 LT=0 


114 64 


115 


72 1 84 


102 


Summer 


LT=1 2 LT=0 


134 77 


150 116 


83 89 


120 M10 


LT=12 I LT=0 


128 66 


138 1123 


Winter 


LT=12 I LT=0 



Interpolation Rules: R (12-months-running mean value): Linear interpolatfon 

LT (local time): Day-night transitions; see section 3.1 .2.3, page 48. 

Modip (modified dip p Equation [3.8]): EPSTEP (p; B 0 (18), B 0 (45), 30, 1 0); see section 

3.1 .2.1 , page 46. B 0 is kept constant for |p| < 18. 


Gulyaeva (1987) has established, based on ionosonde data the following 
relationship between the F2-peak height, hmF2, and the height, ho.5. where 
the density profile has decreased to half the F2-peak density 


P= h a 5 /hmF2 = EPSTEP (2; 0.6, 0.8, 20s, 15) 


(3.17) 


with the seasonal parameter s taking the values 1, 2 3 for winter, equinox, 
and summer, respectively, or using the day number d approximately 

s(d) = 2 - cos (27t/365 • d) (3.18) 


56 










































Fig. 4. Bottomside function for various B 1 values. 


Assuming a bottomside profile as given by Function (3.16), one can use 
Gulyaeva's model to determine the thickness parameter Bo (Bilitza and 
Rawer, 1990). 

B 0 = hmF2 [1 - p)/C (3.19) 

where C is a function of Bj and takes the values 


Br 

3 

3.5 

4 

4.5 

5 

C 

0.755566 

0.778596 

0.797332 

0.812928 

0.82614 


Gulyaeva's (1987) model has also shown good agreement with incoherent 
scatter measurements at Arecibo (Mahajan and Kohli, 1987). 


3.2.3 FI-Layer ( hmFl to HZ) 

Ducharme et al. (1971, 1973) established, based on a large amount of 
ionosonde data, the following relationship for the variation of the FI plasma 
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frequency foFl with solar zenith angle X, solar activity (R 12 ). and magnetic 


dip latitude (T). 

JoFl = Js cos n X (3.20a) 

f s =/o {fwo -So) R 12 / 100 (320b) 

Jo = 4.35 + 0.058 m - 0.00012 T 2 (3.20c) 

f 100 = 5.348 + 0.011 IT I - 0.00023 T 2 (3.20d) 

n = 0.093 + 0.0046 I T I - 0.000054 T 2 + 0.0003 R i2 (3.20e) 


Their model provides also a critical solar zenith angle X s for the occurrence 
probability of the Fl-feature. An Fl-layer is only assumed to exist when the 

solar zenith angle is smaller than X s . 


X s = Xo (*100 - *o) R 12 /IOO 

(3.21a) 

*0 = 49.85 + 0.35 1 T 1 

(3.21b) 

*100 = 38.96 + 0.51 ITI 

(3.21c) 

IRI omits the Fl-feature at night and In winter in accordance with the 
experimental evidence. 


The FI -peak height. hmFl, is found as the height at which the bottomside 
IRI profile reaches the FI -peak density. It Is important to notice that hmFl 
will be affected by the choice of the bottomside thickness parameter Bo. for 
which in IRI-90 two options are given, as described in the previous section. 
Bilitza and Rawer (1990) indicated that option (2), Gulyaeva’s (1987) model, 
produces hmFl values closer to other estimates for this parameter (e.g.. 
Formula [3.331). 

If the FI feature exists, a parabolic Fl-layer is added to the bottomside 
profile 

N [h)/NmF2 = exp (-x^J/coshW + Cj {[hmFl - h)/B 0 ) 1/2 (3.22a) 


with 


fEPSTEP (18; 0.09, 0.2, 30, 10) for In I < 18 ( 

c i = [EPSTEP (|x; 0.09, 0.2, 30, 10) else 

where p is the modified dip latitude (Equation [3.8]) and the EPSTEP 
function is given in Equation (3.2). 


58 


3.2.4 Intermediate Region (HZ to hyi9 

This Is the region where the upper profile (normalized to NmF2 ) is merged 
with the lower profile (normalized to NmE). First, the height HST is found 
where the FI profile function reaches the density value of the E-peak, NmE 
(see Figure 3). It may occur that the profile function does not reach down to 
a very low NmE value or that the height HST is found below the top of the E- 
valley. In such cases the parameter Bj, In the bottomside profile function 
(3.16a) is changed from 3 to 3.5 to 4 and finally 4.5 (see Figure 4). If HST 
cannot be found with these Bj values, the gap between HZ and the E-valley 
top is closed by linear interpolation. The upper boundary of the 
intermediate region Is defined by 

HZ = ( h x + HSTi/2 (3.23a) 


with 


jhmFl if FI -layer Is present 
“ 1 hmF2 else 


(3.23b) 


Starting from HZ the profile is bent downward parabolically such that it 
meets the E-valley top hvr- This is accomplished with the altitude 
transformation 

h’ - HZ = T/2 ± (J2/4 - T (h - HZ)) 1 ' 2 (3.24a) 


where 


T = (HZ - HST)2/(HST - hj) (3.24b) 

and 


hj 


hvr 

hmE 


if valley exists 
else 


(3.24c) 


hvT is the height at the top of the valley and hmE is the height of the E- 
peak. 


3.2.5 E-Peak and Valley (hvr to hmE) 

The E-peak critical frequency is described by the model developed by 
Kouris and Muggleton (1973a, b) for CCIR (1973). It consists of four factors 

fo& = A»B>C»D (3.25a) 
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depending on the 12 -months-running mean of the solar 10.7 cm radio flux 
(COV 12 ), on season (X n oon is the solar zenith angle at noon), on geodetic 
htitudc ftp ). and on the solar zenith angle (X). 


A = 1 + 0.0094 (COV 12 

- 66) 

(3.25b) 

B = COS m ^noon 



(3.25c) 

f-1.93 + 1.92 coscp 

for 

1 cp 1 < 32° 

(3.25d) 

m = to. 11— 0.49 coscp 

for 

1 cp 1 > 32° 


f23 +116 coscp 

for 

1 cp 1 <32° 

(3.25e) 

C = 192 + 35 coscp 

for 

1 cp 1 £ 32° 


D = C0S n 2 a 



(3.25f) 

n -I 1 * 2 

for 

1 cp 1 >12° 

(3.25g) 

n "ll.31 

for 

1 cp 1 £ 12° 



Different from CCIR (1973), a modified zenith angle X a was introduced to 
improve the nighttime variation (Rawer and Bilitza, 1990). 

x a = X - 3 EPS-1 ( X\ 89.98, 3) (3.25h) 

During all times foE is kept above or equal to the observed minimum value 

foEmin = 0.121 + 0.0015 [COV 12 - 66) (3.26) 

A constant peak height hmE = 105 km is assumed for all conditions. 

The F-region peak and valley profile is described by a fifth order power 
series 

N{h)/NmE = 1 + Ej x 2 + E 2 x® + E 3 x 4 + E 4 x 5 (3.27a) 


with 


x = h - hmE 


(3.27b) 


Parameters E lt E 2 . E 3 , and E 4 are determined from the valley width, HER, 
the valley depth, DP , the distance between valley base and hmE, HABR, and 
the logarithmic derivative at the valley top, DLN. The parameter values 
listed in Table 4 were obtained from incoherent scatter measurements. 
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Table 4: Valley Parameters 



HABRi km 

HBRl km 

DP/% 

IHZ2BBB 

hmE/km 

PHal 

Modip =18 

28 

45 

81 

0.06 

105 

Modip = 45 

28 

67 

81 

0.06 

105 

Noon 

Modip = 18 

0 

0 

0 

0 

105 

Modip = 45 

10.5 

17.8 

winter: 10 
else: 5 

summer: 0.01 
else: 0.016 

105 


Notes: DP * 1 00 (NmE - A/yg)/A/mE; N^q is the density at the valley base. 


Interpolation procedure is the same as indicated in Table 3 for B 0 . 

For the deep nighttime valley the power series is fitted to the logarithm of 
the electron density rather than to the density itself. At low and equatorial 
latitudes the valley feature has been observed irregularly during daytime, but 
not frequently enough to be included in a monthly average profile (e.g., 
Mahajan et al., 1990). In the case of zero valley depth (no valley), the FI- 
layer function is merged all the way down to hmE. 

Since it appears quite irregularly, the sporadic E phenomenon is not 
Incorporated in the IRI profile description. 


3.2.6 D-Region and E-Bottomside (hmE to HA) 

Most D-region profiles exhibit a characteristic inflection point with the 
height hmD and the density NmD. Mechtley and Bilitza (1974) found the 
following empirical relationship from rocket measurements (see also Bilitza 
1981a). 

NmD/ ( 1 0 8 m~ 3 ) = (6.05 + 0.088 R 12 ) exp(- 0.1/cos 2 - 7 *) (3.28) 

where R 12 is the 12-months-running mean of the solar sunspot number and 

X is the solar zenith angle. NmD is kept greater or equal to a minimum 
(nighttime) value of 4 • 10 8 m -3 . 

The density profile starts at height HA and is represented by a third order 
polynomial (Mechtley and Bilitza, 1974; Bilitza, 1981a). 

N [h)/NmD = exp (Fj x + F 2 x 2 + F 3 x 3 ) (3.29a) 


with 


x= h - hmD (3.29b) 

Different F 3 parameters are used above (F 3 a ) and below (F 3 b ) the inflection 
point to account for the drastic change in scale height. The parameters 
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Table 5: D-Region Parameters 


Time 

Geographic 

Latitude 

hmD/km 

HA/ km 

F,/km 

FJ km* 

F 3 */km 3 

Fj b /km 3 

hDX/ km 

Day 

Low 

80 

65 

0.02 

-2*1 0- 4 

9.37-1 0 3 

4.89-1 O' 4 

85.6 

Middle 

80 

65 

0.05 

-1.25-1 O' 3 

8.18-10" 3 

1 .707-1 0- 4 

85.6 

Night 

All 

88 

80 

0.05 

-1.25-10" 3 

8.79-1 O’ 3 

1 .22-1 O' 2 

92.5 


listed in Table 5 were obtained from the rocket measurements compiled by 
Mechtley and Bilitza (1974). 


Somewhat arbitrarily an exponential function is used to connect the D-region 
profile with the E-layer peak density NmE and height hmE. 

N [h)/NmE = exp(-D] {hmE - h)*) (3.30) 

The parameters Dj and K are determined such that the function agrees with the 
D-region profile function and its first derivative at the height hDX (see Table 5): 

K = -DN { hmE - hDX)/(NDX In [NDX/NmE)) (3.31a) 

Dj = DN/{NDX • K {hmE - hDX) K" 1 ) (3.31b) 

NDX is the electron density at hDX and DN the derivative dN/dh at that height. 

Comparisons with radio wave propagation data have resulted in several, 
sometimes conflicting, proposals for changes of the IRI D-region profile (Singer 
et al., 1984; Ferguson and McNamara, 1986; Pintado et al., 1987; Oyinloye, 
1988). Since, however, these are all indirectly deduced profiles based on 
certain assumptions about the collision frequencies, additional experimental 
evidence is needed before IRI can be changed with confidence (Serafimov et al., 
1985). 


3.2.7 LAY Functions for Middle Ionosphere (hmE to hmF2) 

A combination of four LAY functions (see section 3. 1.2.2, page 48, Equation [3.5]) 
can be used to represent the electron density profile from the E- to the F2-peak 
in analytical form (Rawer, 1984, 1986, 1988; Suchy and Rawer, 1988). 

4 

log(N( h)/NmF2) = £Aj LAY(h; hmF2, HXj, SCj ) (3.32) 
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A preliminary set of parameters HXj, SCj is listed in Table 6 (Bilitza and 
Rawer, 1990). The amplitudes Aj are obtained automatically by least-square- 
fitting to the constraints listed in Table 7. All peak and valley parameters 


Table 6: Parameters for the Four LAY Functions 


■ 

HX 

SC 

■ 

Day 

Night 


1 

0.9 hmF2 

^0.5* 

0.1204 ( hmF2-h 05 ) + 56.8 

2 

hmFV 

(hmF2+h VJ )/2 

10 

3 

hvB 

hvB 

9 

4 

hmE- 6 

hmE 

6 


* 

(hmF2 + h 0 5 )/2, if FI -layer is not present, 
t h 05 \s the height, where N(h) = 0.5 NmF2{see Equation [3.17]). 


Table 7: Constraints for LAY-Formalism 




h 

N 

dN/dh 

Welaht 

Day 

1 


Night 

1 

i 

ho.5 

0.5 NmF2 

- 

fit 

1 

2 

2 

h VT 

NmE 

- 

1 

1 

3 

3 

hvB 

Nvb 

- 

2 

3 

4 

4 

hmE 

NmE 

- 

5 

5 

5 

5 

^BE* 

Nvb 

- 

1 

0.5 

6 

- 

hmFI 

NmFI 

- 

3t 

- 

7 

6 

hvB 

- 

0 

50 

50 

8 

7 

hmE 

- 

0 

500 

500 


* 

hgE = hmE— {hyg— hmE). 
i Zero, if FI -layer is not present. 

tt f = EPS 0 (z; 0.15, 0.1 ) with z - log (N 0 jNmF1) and N o s - 0.5 NmF2. 

= 1 if FI -layer is not present. 

= 0-5 if (NmFI - N 0 S )(hmF1 - /) 05 ) < 0. 
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are the same as those described for the standard IRI-90 profile, except for 
hmFl which is calculated with the widely used formu a 

hmFl /km = 165 + 0.6428 X (3.33 


X is the solar zenith angle in degrees. Because there are more constraints 
(seven or eight) than free parameters (four), the degree to which the 
individual constraint is enforced will depend on the weight that is assigned 
to it during the fitting procedure. For more information and applicabo t 
other parts of the density profile, see Rawer and Bilitza (1990). 


3.2.8 Plasmaspheric Extension 

Rvcroft and Jones (1985, 1987) have developed a plasmaspheric extension 
for the IRI electron density model. Their diffusive equilibrium model 
describes the field-aligned distribution of plasma out to L-values oi 6 A 
merging procedure with the IRI topside profile is suggested at a reference 
level near 650 km. This plasmaspheric extension is not yet included in IRI 
90 but may be available from the authors. 


3.2.9 Ionospheric Electron Content 

McNamara (1983, 1984, 1985) has compared a large data base of total 
electron content (TEC) observations with the TEC values computed with 
IRI-79 He finds generally good agreement at midlatitudes; however. IRI-79 
underestimated the TEC observed close to the magnetic equator by up to a 
factor of 2 The improved topside profile in IRI-86 helped to reduce this 
discrepancy to about 10% to 20%. Using the new option for the bottomside 
thickness parameter Bq will further improve the agreement. 


Application of IRI 
al. (1988) and by 
measurements is 


for satellite orbit determination was discussed by Bilitza et 
Coster et al. (1990). A recent comparison of IRI and TEC 
given in section 2.4, page 33, of this report. 
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3.3 Plasma Temperatures 


The IRI-90 profiles for the electron and ion temperatures extend from 120 km to 
1000 km. First order estimates can be obtained up to 3000 km; above 1000 km 
the two temperatures are assumed to coincide. 

At the lower boundary, HTA =120 km, thermal equilibrium is assumed between 
neutrals, ions, and electrons 

T e(HTA) = Ti (HTA) = T n (HTA) (3.34) 

in accordance with most of the observations. Differences between electron and 
neutral temperatures at these low altitudes have been reported by Oyama et al. 
(1980). 

With increasing altitude the electron temperature is always kept above (or equal 
to) the ion temperature, which in turn is kept above the neutral temperature 

T e (h) > T,(h) > T n (h) (3.35) 

Variations with solar activity are not explicitly included in the IRI-90 plasma 
temperature models. The above stated constraints (coupling to neutral 
temperature) can, however, cause solar activity variations when all temperatures 
are close together. In such cases, for example, at nighttime, the plasma 
temperatures vary with solar activity in the same way as the neutral temperature 
does. A review of observed and modeled solar activity variations of plasma 
temperatures was given by Bilitza and Hoegy (1990). 

IRI uses the COSPAR International Reference Atmosphere (CIRA) as its neutral 
temperature model. The IRI computer code contains subroutines to calculate the 
CIRA (1972) temperatures. It is planned to incorporate the newly released CIRA 
(1986) model into IRI, as soon as an appropriately simplified version is available. 


3.3.1 Electron Temperature 

The development and improvement of the IRI electron temperature model 
was described by Bilitza (1981b, 1985b) and Bilitza et al. (1985). The electron 
temperature profile is subdivided into six regions defined by the region 
boundaries AHHj [j = 0 to 6) as shown in Figure 5. Using Booker's (1977) 
approach, the temperature profile is approximated by a straight line in each 
region. Thus, the temperature can be described as (see section 3. 1.2.1, page 
46, Equation [3.3]) 

T e (h) = ATE 0 + B(h; AHH 0 , 5, ST, AHH, DTE) (3.36a) 

with the region gradients 


STj = (dT e /dh)j_i, j = ( ATEj - AT^_i)/(AHHj - AHHj_i) (3.36b) 
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Fig. 5. Buildup of the IRI electron temperature model. 


The region boundaries AHHj, the electron temperatures ATEj at these 
boundaries, and the transition thicknesses DTEj are given in Table 8. If we 


Table 8: Electron Temperature Parameters 


i 

AHHji km 

DTE/km 

ATE/ 

Source Model 

Data Base 

*12* 

0 

120 

- 

T n (120km) 

CIRA (1986) 

Incoherent 
scatter, satellites 

T n varies 
with R 12 

i 

fom 

5 

TE m 

- 

Incoherent scatter 

60-120 

m 

300 

5 

7£ S7 <300km) 

Brace and 

AE-C 

10-20 

1 

400 

10 

7£ ST <400km) 

Theis (1981) 

AE-C 

20-50 

H 

600 

20 

7£ SP (600km) 

Spenner and 
Plugge (1979) 

AEROS 

20-40 

Bl 

1400 

20 

7£ S rt 1400km) 

Brace and 

ISIS 2 

40-70 

H 

3000 


7£ n r(3000km) 

Theis (1981) 

ISIS 1 

70-110 


* R 12 is 12-months-running mean solar sunspot number. 
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(3.37) 


evaluate function (3.36) at the altitudes AHHj, we find 

T e (AHHj) - ATEj = DSSj » ( ST J+1 - STj) DTEj ln(2) 

DSSj describes the difference between the skeleton profile and the 
smoothed function at AHHj. The ATEj values are, however, well defined 
temperatures, which should be on the final profile. This can be 
accomplished to a good degree by introducing new temperatures 

ATE/ 1 ) = ATEj - DSSj (3.38) 

and using these as the anchor points for the skeleton profile. 

Table 8 lists also the source models used to obtain the temperature values 
ATEj and the data base behind these models. Brace and Theis (1981) apply 
spherical harmonics analysis to describe their AE-C and ISIS satellite data in 
terms of magnetic dip latitude (Equation [3.11]) and local time. Sets of 81 
coefficients are provided for equinox and solstice at 300 km, 400 km, 1400 
km, and 3000 km. The model by Spenner and Plugge (1979) is based on 
AEROS satellite data and describes variations with geomagnetic latitude, 
longitude, and altitude. Two sets of coefficients are given, one for 3 a.m. 
local time and one for 3 p.m. (both for winter). Omitting the small 
longitudinal variation, Bilitza (1981b) obtained the following approximate 
expressions. 


TEsp (600 km, 4>)/K = 
(3.39) 


2900-5600 EPSi (4>; 0, 11.35) 
EPSTEP (I I ; 839, 2000, 45, 5) 


at 3 p.m.l 
at 3 a.m. J 


4> is the geomagnetic latitude in degrees. Local time variations are 
described with the day-night transition function explained in section 
3. 1.2.3, page 48 (Equation [3.6]). 

During daytime at low latitudes, the temperature profile frequently exhibits 
a local maximum at about 270 km. The height h m and temperature TE m of 
this feature were obtained from incoherent scatter measurements. 

h m /km = 210 + 60 exp (- (4>/22.4) 2 ) (3.40a) 

TE m /K = 1500 + 800 exp (- ($/33) 2 ) (3.40b) 

During nighttime TE m is set equal to the neutral temperature at h m = 150 km. 

The IRI model has been compared with a large number of rocket, satellite, 
and incoherent scatter radar measurements. In general good agreement 
was found; discrepancies are mostly due to the Insufficient description of 
changes with season and solar activity (Bilitza and Hoegy, 1990; Buonsanto, 
1989). 
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IRI offers a second option for the determination of the temperatures 2 . 
ATE? and ATE4. They can be calculated with a model that makes use of the 
strong anti-correlation between electron temperature and density in the 
daX® c tonosphere (Brace and Theis. 1978; Bllitza et ah. 1985) In this 
case” users have to provide the electron density at the heights at which they 
choose to use this option (AHH,, AHH 2 . and/or AHH 3 ). The model by Brace 
and Theis (1978) describes the electron temperature in terms of electron 
density and altitude. It was extended by Bilitza et al. (1985) to include solar 
activity variations. 


3.3.2 Ion Temperature 

Similar to the electron temperature, the ion temperature profile is 
represented by a Booker-function 

T t (h) = ATIo + B(h; XSM 0 , 3, STI, XSM, DTI) 


with the gradients 

STIj = [ATIj - ATIj-i) / [XSMj - XSMj-j) (3 41b) 

the section boundaries XSMj, the ion temperatures at these boundaries ATIj, 
and the transition thicknesses DTIj (see Table 9). 


Table 9: Ion Temperature Parameters 


i 

XSAf/Vkm 

ATli 

DTIj 

0 

HS 

Tn(HS) 


1 

430 

TIM 

10 

2 

HTT 

T e (HT7) 

10 

3 

AHH 5 ‘ 

ATEs* 

20 

4 

AHH a* 

ate 6 * 

- 


If HTT> AHH 5 , then XSM 3 = AHH 6 and ATI 3 - ATE 6 . (AHH, and ATE , are electron temperature parameters defined 
in section 3.3.1 , page 65. and Table 8.) 


Incoherent scatter observations of ion temperature show a distinct change 
in gradient at about 430 km. Using AEROS satellite measurements the 
latitudinal variation of ion temperature at this point has been described as 


(Bilitza, 1981b) 
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TIM = 


(3.42a) 


1240- 1400 EPSi (*; 0. 11.11) day 1 
1200 - 300 cos^Z night] 

with the geomagnetic latitude 4> (in degrees) and 

Z = 0.47 I <t> I + 0.024 4> 2 (3.42b) 

The height HS is found as the point on the neutral temperature profile, 
whose tangent meets the point (XSMj/T7M). Below HS, T, = T n is assumed. 

HTT is the height at which the ion temperature profile meets the electron 
temperature profile, assuming a constant ion temperature gradient STI 2 
above XSMi. STI 2 is chosen as 

ST '* = {0 5^ht} (3.43) 


based on incoherent scatter data. Above HTT, Tj = T e is assumed. 

Comparisons with incoherent scatter data and with satellite measurements 
have indicated the need for a closer coupling between ion and neutral 
temperature at intermediate altitudes (Bilitza and Hoegy, 1990). See also 
Buonsanto (1989). 
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3.4 Ion Composition 

IRI describes the relative percentage ion densities (ion composition) R lon (h); 
thus, the absolute ion densities are given by 

Nion(h) = N e(h) Rion(h) (3.44) 

The model describes the percentage of 0 + and 0 2 + ions and fills up to 100 % 
with light ions in the topside and with NO + ions at lower altitudes. A tenth 
of the light ions is assumed to be He + , the rest H + . 

Booker's (1977) approach is used to represent the 0 + and 0 2 + profiles (see 
section 3. 1.2.1, page 46). For 0 + ions the profile function reads (see 
Equation [3.3]) 

RO(h) = RO 0 • exp( B(h; HO 0 , 4, MO, HO, DO) ) ( 3 . 45 ) 

with the thicknesses DOj = 9, 5, 5, 50 km for j = 1 , 2, 3, 4. At low altitudes 
the gradients MOi, M0 2 and the transition heights HOj, H0 2 were obtained 
from the compilation of rocket data (ion mass spectrometer) by Danilov and 
Semenov (1978). The variation of these parameters with solar zenith angle 
X was approximated by the step function (3.2) 

MOi, MO 2 , HOj = EPSTEP (z; A, B, z s , dj (3.46) 

where z = cos X. Because of the limitations of the data base the model is 
restricted to sunlit hours. In IRI, z is set equal to 0 for X greater than 90°. 
Parameters A, B, z s , and d z vary with season and solar activity as shown in 
Table 10. The topside part of the 0+ profile is yet relying on a very small 
data base of satellite data (AEROS, OGO 6 ; retarding potential analyzer). It 
distinguishes only day and night conditions and two latitude classes as listed 
In Table 1 1 . 

Only two Epstein functions are needed to represent the 0 2 + profile. 

R02(h) = RO2 0 • exp( B(h; HO2 0 , 2, M02, H02 , D02)) (3.47) 

with the transition thicknesses D 02 j = 5 km and D02 2 = 5 km and the 
gradients and transition heights as listed in Table 12. 

The following adjustments are made automatically during the program 
execution: 

( 1 ) The height HO 3 is found where the upper 0 + skeleton profile (defined 
by the point (HO s /l%), the transition height HO 4 , and the gradients 
MO 4 , MO 5 ) reaches 100%. 

HO 3 = H0 4 + ( ln(100) - MO s ( HO 4 - HO 5 D/MO 4 (3.48) 


71 

PRECEDING PAGE BLANK NOT FILMED 



Table 10: Parameters for the Lower Part of the IRI 0 + Model 



Winter 

P Summer : 

Equ 

Inox 

cov < 100 

COV > 100 

mu 

cov > 100 


A 

0.13027 

0.161 

0.092 

mm 

1 

AfO, 

B 

0.08018 

0.069 

0.072 


Krai 

%/km 

Zs 

0.25 

0.18 

0.21 

SP mi 

PH 


d z 

0.04216 

0.254 

0.014 

H39H 



A 

0.00999 

0.0216 

0.03863 

0.0195 

0.01 

M0 2 

B 

- 0.00686 

0.0161 

0.01389 

0.0067 

0.00127 

%/km 

z s 

0.1 

0.1 

0.12 

0.1 

0.09 


dz 

5.113 

0.03014 

0.05762 

0.04 

0.05 


A 

180 

167 

168 

172 

185 

HO, 

B 

170 

152 

165 

158 

167 

f 

km 


0.15 

0.17 

0.258 

0.24 

0.18 



0.1175 

0.04916 

0.008 

0.01 

0.015 

HOzI km 


290 

290 

237 

290 | 

HOolkm 

(H0 2 +H0 3 )l 2 

RO 0 /% 

| 100 


Table 11: Parameters for the Topside Part of the IRI O* Model 



1 !Latitudel<30 1 

|Latitude|>30 1 


Day 

Niqht 

Day 

Night 

HCVkm 

695 

570 

695 

575 

HO s *l km 

2177 

1040 

3367 

1380 

m&BX5SM 

• 

M0 4 l%km"1 

-0.000781 

-0.002 

-0.000786 

-0.00126 


0.00264 

-0.0052 

-0.00165 

-0.00524 | 


* HO s is the height at which the percentage of 0 + ions has decreased to 1%. 


HO 3 has to be above H0 2 (see Table 10), which is the height at which 
the lower 0 + skeleton profile reaches 100%. If HO 3 is below HO 2 , then 
the gradient MO 4 is successively decreased by 0.001 km' 1 until HO 3 is 
above H0 2 . 

(2) Using the reference point H02ol / R02o and the gradients M02 from 
Table 12, a new reference point is determined at the height ( H02 2 + 
H02 3 )/ 2 and used in Formula (3.47), thus making sure that the 
reference (anchor) point for the 0 2 + profile is not located in one of the 
transition regions (see section 3. 1.2.1, page 46). 
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Table 12: Parameters for IRI Ol Model 




Winter 

Summer 

l_ Equinox | 



cov< 100 

COV> 100 

cov< 100 

COV> 100 

WHEM 

a 

136 

130 

130 

125 

125 

H *E 

b 

-9 

-10 

-10 

0 

0 

IMAiUm 

a 

181 

198 

190 

182 

170 

II km 

b 

-26 

0 

0 

0 

0 

\MZnHM 

a 

0.02994 

0.05922 

0.05107 

0.0307 

0.02806 

1 %/km 

b 

-0.04879 

-0.07983 

-0.07964 j 

-0.04968 

-0.04716 


a 

-0.01396 

-0.00397 

0.00097 H -0.00248 I 

-0.00066 

II %/km 

b 

0.00089 

0.00085 

-0.01118 

-0.02451 

-0.02763 


a 

-0.09929 

-0.00313 

-0.02614 

-0.00313 

-0.02247 

II %/km 

b 

0.05589 

0 

-0.09537 

0 

-0.01919 


140 

160 

140 

140 | 

140 

|| R02o/% | 

45 

31 

30 

37 

37 | 


. . , H02\ , f cos X . . 

No,e: M02 >-a + bz z ” \o f or * > 9 0 * ls solar zenith an 9 le - 

(3) The height HO max is determined where the 0 + profile function (3.45) 
reaches its maximum ROmax- If ROmax is greater than 100%, then RO 0 
is reduced successively by 0.01%, until PO^^is less or equal to 100%. 
If at this altitude the percentage of O 2 + Ions (Function [3.47]) is larger 
than 100 - ROmax* then the gradient MO 23 is decreased successively by 
0.02 km -1 , until the sum of 0 + and C> 2 + percentages is less than 100%. 

(4) At the height HO max the ratio 

r = (100 - ROmax ~ R02{H0 ma J)/R02{H0 ma J (3.49) 

is determined and is used to calculate the percentage of NO + above 
HO max . Herewith, interference between the molecular ion profiles and 
the light ion profiles is avoided. 

Finally, the percentages of NO+, H + , and He + ions are determined: 


|l00 - RO(h) - R02(h) hZHOnax 
RNOln) - j r R02(h) h > H<W 


(3.50) 

f(100 - RO(h) - R02(h) - RNO(h))/(l - P) 
RH(h) = | 0 

h> HOmax 
h < HOmax 

(3.51) 

RHE(h) = RH(h) P/(l - P) 


(3.52) 
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P = 0. 1 is the ratio of He + ions to the sum of all light ions (H + and He + ). 

Comparisons with incoherent scatter data (Bilitza, 1984) with OGO 6 data 
(Kutiev et al., 1984a, b) and with ISS-b data (Goel and Rao, 1984) showed 
that IRI overestimates the percentage of light ions in the topside 
ionosphere. Discrepancies were also found at low altitudes in comparison 
with S3-1 satellite data (Philbrick and Bhavnani, 1983; Philbrick et al., 

1984) and in comparisons with equatorial rocket data (Sridharan et al., 

1985) . A review of empirical ion composition modeling and of possible 
improvements of IRI was given by Bilitza (1990). 

An alternative model was proposed by Danilov and Yaichnikov (1985). It 
describes the relative densities of 0 + , H + , He + , NO + , 02 + , N + , and Cluster ions 
in the altitude range from 75 km to 1000 km, taking into account variation 
with solar zenith angle, latitude, season, and solar activity. It is based on a 
compilation of satellite (Electron 2, 4, S3-1) and (high apogee) rocket data. 
IRI-90 offers this model as a second choice for the computation of the 
ionospheric ion composition. A listing of the computer code follows below. 
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subroutine ioncom(h , z , f , f s , t , cn) 

c 

c ion composition model 

c A.D. Danilov and A.P. Yaichnikov, A New Model of the Ion 
c Composition at 75 to 1000 km for IRI, Adv. Space Res. 5 #7 

c 75-79, 107-108, 1985 
c 

c h 

c z 

c f 

c f s 

c t 

c cn(l) 

c cn(2) 

c cn(3) 

c cn(4) 

c cn(5) 

c cn(6) 

c cn(7) 


dimension cn(7) ,cm(7) ,hm(7) ,alh(7) ,all (7) , beth (7) , 

betl (7) , p(5,6,7) ,var(6) ,po(5,6) ,ph(5,6) , 
pn (5,6), phe (5,6), pno (5,6), po2 (5,6), pci (5 , 6) 


altitude in km 

solar zenith angle in radians 
latitude in radians 
10.7cm solar radio flux 
season (month) 

0+ relative density in percent 
H+ relative density in percent 
N+ relative density in percent 
He+ relative density in percent 
N0+ relative density in percent 
02+ relative density in percent 
cluster ions relative density in percent 


k 

k 

k 

k 


k 


k 

k 

k 


data po/ 4*0. ,98.5,4*0. ,320. ,4*0. ,-2.59E-4,2.79E-4,-3.33E-3, 
-3 . 52E-3 , -5 . 16E-3 , -2 . 47E-2 , 4*0 . , -2 . 5E-6, 1 . 04E-3 , 

-1 . 79E-4 , -4 . 29E-5 , 1 . 01E-5 , -1 . 27E-3/ 
data ph/-4. 97E-7 , -1 . 21E-1 , -1 . 31E-1 ,0. ,98. 1 , 355 . , -191 . , 

-127. ,0. ,2040. , 4*0 . , -4 . 79E-6 , -2 . E-4 , 5 . 67E-4 
2 . 6E-4 , 0 . , -5 . 08E-3 , 10*0 . / 

data pn/7 . 6E-1 ,-5.62,-4.99,0. ,5.79,83. ,-369. ,-324. ,0. ,593. , 
4*0. , -6 . 3E-5, -6 . 74E-3, -7.93E-3, -4. 65E-3,0. , -3.26E-3, 
4*0. ,-l . 17E-5 , 4 . 88E-3, -1 .31E-3, -7.03E-4,0. , -2. 38E-3/ 
data phe/-8 . 95E-1 ,6.1,5.39,0. ,8.01,4*0. ,1200. ,4*0. , -1 .04E-5 
1 . 9E-3, 9 . 53E-4 , 1 . 06E-3,0. ,-3.44E-3,10*0./ 
data pno/-22. 4, 17. 7, -13. 4, -4. 88, 62. 3, 32. 7,0. ,19.8,2.07,115. , 
5*0. ,3. 94E-3,0. , 2 ,48E-3,2. 15E-4, 6. 67E-3, 5*0. , 

-8 . 4E-3 , 0 . , -3 . 64E-3 , 2 . E-3 , -2 . 59E-2/ 
data po2/8. ,-12.2,9.9,5.8,53.4,-25.2,0. ,-28.5,-6.72,120. , 
5*0. , -1 . 4E-2 , 0 . , -9 . 3E-3,3 .3E-3.2. 8E-2,5*0. ,4.25E-3, 
0. ,-6.04E-3,3.85E-3,-3.64E-2/ 
data pcl/4*0 . ,100. ,4*0. ,75. ,10*0. ,4*0. ,-9.04E-3,-7.28E-3, 
2*0. ,3.46E-3,-2. 11E-2/ 


DO 8 1=1,5 
DO 8 J=1 , 6 

P(i, j ,l)=po(i, j) 
p(i,j ,2)=ph(i, j) 
p(i. j , 3)=pn (i , j ) 
P(i. j ,4)=phe(i, j) 
P(i> j ,5)=pno(i, j) 
P(i, j , 6) =po2 (i , j ) 
P(i, j ,7)=pcl(i, j) 

8 continue 


s=0 . 
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7 


1 

2 

3 

4 


5 

6 


c 

c 


do 5 i=l ,7 

do 7 j=l ,6 

var(j) 

& 

& 


p(l, j,i)*cos(z) + p(2, j ,i)*cos(f) + 
p (3 , j , i) *cos (0 .013* (300. -f s) ) + 
p (4 , j , i) *cos (0.52*(t-6.) ) + p(5,j,i) 


continue 
cm (i) = var(l) 

hm (i) = var (2) 

all (i) = var (3) 
betl (i) = var (4) 
alh(i) = var(5) 
beth(i)= var(6) 
hx=h-hm (i) 
if(hx) 1,2,3 

cn (i) = cm(i) * exp( hx * (hx * 
goto 4 

cn (i) = cm(i) 
goto 4 

cn (i) = cm (i) * exp( hx * (hx * 
continue 

if (cn (i) . LT. 0.005* cm (i)) cn(i)=0. 
if (cn (i) . GT . cm(i) ) cn(i)=cm(i) 
s=s+cn(i) 
continue 
do 6 i=l,7 

cn(i)=cn(i) /s*l00. 


all (i) + betl (i) ) ) 


alh(i) * beth(i)) ) 


return 

end 
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3.5 Ion Drift 

A model of the ionospheric ion drift was developed for IRI-90 by 
Kazimirovsky and his colleagues based on ground-based drift measurements. 
It provides the horizontal E- and F-region drifts (zonal and meridional) for 
different geomagnetic latitudes, local times, seasons, and solar activities. 

The model is described in section 2.5, page 41, of this report. A listing of 
the computer code and coefficient files follows below. The program and 
coefficients may be available in computer-accessible form from the model 
authors. 
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A Computer Code for Determination VX, VY 


PROGRAM GlD3(lNPUT,0jT PUT, FARl,FAR2>; 

£***„****„**#***********+*****************************************) 
(* TEST CALCULATION OF IONOSPHERIC DRIFTS GLOBAL M 0 0 E L * > 

(************************************************************** ,t **^ 

TYPE T 1 = ARRAY (.1,,15,) OF REAL? 

T 2 = ARRAY (.5. .10,0. ,23,) OF REAL; 

VAR F A R 1 , F A R 2 : TEXT,' 

3F AR : T 1 j 

0 X , 0 Y : T 2 ; 

I,J,PMAX,REG.P : integer; 

F I , L T , V < , v Y : real; 

AR : ARRAY (,!., 78,1, .15.) OF R E A |_ ; 


procedure drir«< r c g , p : integer; FI, ITT: real; VAR V X , V Y : REAL >; 

(**.*************#************★***********♦************************) 

(* E , S . KftZ I MI ROVSK Y , 6 , I . ZHOVT Y A N U w . A. CHERNIGOVSKAYA *) 

< * SIBERIAN IN5TITJTE OF TERRESTRIAL MAGNETISM, IONOSPHERE AND *) 
(* RADIO WAvE PROPAGATION. SIBERIAN pEHAPTMENT OF THE U.S.S.R, *) 

< * ACADEMY OF SCIENCES, S 6 !, 0 3 3 IRKUTSK 33 , P.O.BOX 4026 ,U.S,S,R,*> 

I*****************************************************************^ 

*) 


PROCEDURE CALCULATES ZONAL VX AND MERIDIONAL VY 
COMPONENTS OF HORIZONTAL VELOCITY OF THE IONOSPHERIC 
IRREGULARITIES AT THE HEIGHTS OF f - AND F-REGION FOR PERIOD 
NUMBER 3 , GEOMAGNETIC LATITUDE FlcptGP) AND LOCAL MEAN 
TIME LTT(HOJR), IN M/SEC. 


TO CALCJLATE NX AND VY 

FOR E-’REGI ON IT IS NECESSARY Tu TAKE R 

FORF-REGION r 

FOR THE WHOLE PFRIOD OF OBSERVATIONS 

(ONLY FOR F-REGTON) 
FOR W 1 nTER<F 3 -R THE 4 H 0 L E P E R I 0 l» OF 0 8 S E R V 

FOR S 3 RlNG(FOR THF WHOLE PERIOD OF 3 8 S E R V 

FOR S JM w ER<FOR THE WHOLE PERIOD OF 0 b S t R V 

FOR A J T U M N ( F ) * THE WHOLE PERIOD OF OBSERV 


R E G= 1 
R E G = 2 


(* FOR WInT£R<FOR the WHOLE PERIOD OF OBSERVATIONS) P=2 *> 

(* FOR S 3 RlNG(FOR THF WHOLE PERIOD OF OBSERVATIONS) P=3 *) 

< * FOR SJM^ERCFOR THF. WHOLE PERIOD OF OBSERVATIONS) P = A *) 

(* POR A J T U M N ( F ) * THE WHOLE PERlOu OF OBSERVATIONS) P = 5 *) 

(* FOR W I N T £ R ( F 10.7 * ?00 ) P = 6 *) 

(* FOR S 3 RInG( F 10.7 * 200 ) Ps? ° 

(* FOR S J M M £ R ( F 10.7 = ZOO ) p = 8 ° 

(* FOR AJTUMNC F 10.7 = ZOO > p = <? 

(* FOR W I N T E R < F 10.7 = 70 > Pa1 ° *> 

(* FOR S 3 R I N G ( F 10.7 = 70 > Ps11 

(* FOR SJM^ERC F 10.7 8 70 > P=12 *> 

(* FOR AJTUMNC F 10.7 i 70 > P=13 *) 

* ) 

( * 

(* IF THERE IS NO OATA ThEN V X = V Y s p 9 9 ? . 9 *> 

(* ° 
(* ; X > 0 To EAST, VY > 0 TO NORTH *> 

(********»******K-**************** , ******************************* , * ) 
(* SEE "ADV. SPACE RES," V0L,S*#7,95,1C9M9dS)» V 0 L, 8 »*A, 239(1 988)*) 

(* FOR DESCRIPTION OF The ORIGINAL WORK, *> 

( *******************'**************' , ************************** <, **** > 

<* THIS VERSION OF MODEL IS WRlTiFN IN PASCAL 8000 * > 

<* SEE "PASCAL B 0 0 0 . REFERENCE M A AT J A L (VERSION 2 , 0 ) " , *) 

(i*************************************************************** 0 
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BEGIN <* DRIRR *) 

RA:rPi/ia; raojs^i/isd; 

(* CHECKING T 4 E CONDITIONS OF INFORMATION AVAILABILITY *) 

L T : = L r r ; 

T3 : I F L T > 2 A THEN 3 E 3 I N lT:«LT-2A; GOTO 1 J END? 

if ( f i > 6 o ) or ( f I < 1 o ) then goto i i ; 
k: = (p,1)*s*’Trjnc(Fi/io); 

Pl:=TRUNC(AR(,< f 1S,)>; 

CASE REG OF 

1: IF (PI sO) OR <(P 1 = 1 ) And ( ( L T < 6 ) ok ( L T > 1 8 ) ) > THEN GOTO 11 J 
2: IF <P1*0> OR <<Pl=1) AND ( ( L T > 6 > AND (LT<13>>> THEN GOTO 11 
END! < * CASE RES * ) 

IP TRJNC(PI/10)sTR'JNC((FI- 1)/10 + 1) T H t N 
3E 51 N ( * 1 O 

FOR I 1 : = 1 TO ? DO 
H E G I N 

I ? : = < 1 1 - i > * ? ; 

vo:=ar(.<#I2-m,>; vi: = ar<,k,i**?.>; v 2 : = a r < , k , i 2 3 , > ; 
V3:=ar(.<»i?>4,>; f1:*ar<,k,i^*5.>; f2:=ar(,<,i2+6.>; 

f ? : = a 3 ( . < , I 2 + 7 , ) ; 

CASE II OF 

1 : v < : s C ( L T , V O . V 1 , V 2 , V J . F 1 f F 2 , f 3 > ; 

2: VY:sC(LT,V0»Vl,V2,V5,Fl,F2fF3) 

END ( * C A S F II * > 

END ( * II * ) 

END ( * 1 * ) 

ELSE 

3EG I N < * 2 * > 

P11 :sTRUNC(4RC,K*1 f l5,)>; 

CASE REG OF 

1: IF ( P 1 1 S 5 ) or ( < P 1 1 = 1 ) and U[r< 6 ) OR < L T > 1 8 > ) > 

THEN GOTD 11; 

2: IF ( P 1 1 = 0 > OP < < P 1 1 « 1 ) AND t ( L T > 6 ) AND ( L T < 1 8 ) ) ) 

THEN GOTD 11 

end; < * case Reg * ) 

C* FILLING OF THE ARRAV F *1 
FDR I 1 : = 1 TD l DO 
BEGIN 


I ? : = 

( 1 1 

- 

1 > * 

?; 













F 1R 

I 3 : 

s 

1 T 

0 

6 

DO 











BE 

G I N 

















J : * 

< 

3 - 1 

) * 

6 

*13; 












vo : 

7 

A R ( 

, J 

t 

12*1, 

) ; 

V 1 

: = ap ( 

* J 9 

I 2 + 2 . 

> ; 

V 2 ; 

= AR ( 

t ^ f 

12 * 3 ,); 


v 3 : 

r 

A 3 ( 

. J 

r 

I 2 * A . 

) ; 

F 1 

: = ar ( 

. j, 

12*5 . 

> j 

F 2 : 

= A R ( 

* J f 

1 2*6 , > ; 


F 3 : 

= 

A R ( 

• J 

f 

12*?. 

) ; 

f ( 

, 13^1 

* > : 

s C < L T 

f vo 

, VI 

# V2, 

^ 5 f 

F 1 , F 2 , F 3 ) 

end; 

( 

* I 

V 

* 

> 












(* the latitudinal interpolation of velocity *> 

CASE II OF 

11 v<:sspf(pi/18,fi*pad>; 

2 ! V Y : a s p F ( H I / 1 3 , F I * » A D ) 

END < * CASF II *) 

E NO 

END; < * 2 * ) 

GOTO 12 ; 

11:VX:=9P99,R; vY;=9999.9; 

12 : 

END! ( * DRIRR * > 
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label 11,12.13; 

CONST P I = 3 . 1 M5926; 

var <# pi , 1 1 , 12, 1 5, 1 * J , ‘’ii * integer; 

LT,RA,RAD.V0.V1rV2,V3. F1,F2#F3; « F A L 

f: ARRAY ( . 0 . , 5 . ) OF REAL? 


* ) 
* > 


FUNCTION SPF< H*, C REAL >: REAL* 

<* FUNCTION CALCULATES SALINE INTERPOLATION OF 
(* lautjdinal variation 
L A 3 E L 1 ; 

VAR G.3.R.Z ; ARRAY (.0,.S.) OF REAL,' 

r , < s i . n u : array ( . D . . 6 « > oe real; 

A , B : ARRAY ( , 0 . , 4 . > OF REAL# 

1,10 : integer; 
x o , < i : real; 

BFjIN < * SPF * > 

5 ( . 3 , > : * 0 #’ R < . 5 . ) : O ; 3 ( , 5 . > : = 1 ; «(.0.):=1*' 

Z<.0.):*3*<F<.1.)-F<,0»>>/HH; 

Z(.5.):*3*(FC,S,>“F<«4»))/R^» 

FDR I : = 1 TO A 0 0' 

BEGIN 

G ( . I , > : * 5 . 5 ; R < , I . ) I = 0 . s ; 

Z < . I . > : = 3 * < R c . I , > * < F C . I ♦ 1 , > - F Y . I , > > / H X ♦ 

3(.I.)*(F(.I.)-F( f I-1.)JYHXJ 

e no ; 

FOR I : = 0 To 5 00 9 < 1 I • > ! *-2 ; 

< s i ( . o . > : = o ; n j < , o . > i s o ; 

FOR I : S o To 5 00 

BEGIN 

nj(.I+1.);*<NJ(.I.>*G(.I.)-2<.I.>)/ ( Q^I' > ” <SI< « I « > * GI ’ I ’ >) 

e no; 

FOR^iUs DOWNTO 0 OO ^( . I . ) ! =<SI l ,P1 ,)*•?( . 1*1 , ) + NU( , 1*1 , ) ! 

FOR I:=0 TO '* 00 

A(.I.>:=-2*IF<.I*1.>“F<.I,)>/HX'* - CR<.1.>* v T<iI + T» >) » 

B(. I . > :=-A(.I,> + CFC.I + 1 ,>-F<.I.»>/HX-v,<,I.> 

E n o ; 

FOR I : = o to 5 00 
3 e G I N 

IF (X>s<n»RAO+HX*I>) A N 0 CX<(10*RA0+HX*U + 1>>) THEN 

BEGIN xi:=x; iO:si; GOTO n end*' 

end; 

l4 SPF:=F(/lo!>«-B(.IO.)*XOt3<.lO.)/MX**3**^ , * , A(,IO,>/IHX*HX)*XO**3 

F N o ; (* S 3 F +) 

FUNCTION C < J J * C 0 . C 1 * C 2 . C 3 , F I 1 , F I 2 , F I 3 ; REAL) : R E A L » 

<* FUNCTION CALCULATES TrtF. VALUE OF VELOCITY FOR *) 

(* THE LATITUDINAL NET ROINTS * 5 

BE c!=C0^C1*SIN(RA*JJ+Fn*RAD>*C2*SlN(2*PA*JJ^FI2*RAD)* 

C3*SIN<3*RA*JJ*FI3*RA0) 

E N 0 ; ( * c * > 
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PROC E D JR E NAMECREG, 3 : INTEGER); 

<* PsnCEOJRF ?)hTS THE NME OF THt PtRIOD FOR *> 

(* DIFFERENT R E 3 And P *> 

VAR II: integer; 

BEGIN 

CASE REG OF 

1 : 4 R I T E ( ' E REGION ' > ; 

2 : WRITEC' F REGION * ) 

E N 0 *’ < * CASE R E 3 * > 


I F 






< REG = 2> 

4 N p 

<R = 1 

) 

THEN 

I I : 

= o ; 

IF 

( R E G = 1 > 

An D 

( 3 = 1 

) 

DR 

< REG=2> 

AlNp 

IR = 2 

) 

THEN 

I I : 

al ; 

IF 

( R E G = 1 ) 

and 

(5-2 

) 

DR 

< ft 6 G = 2 ) 

A N p 

Uc3 

) 

THEN 

1 1 ; 

*2 ; 

I F 

( R E G = 1 > 

and 

<3 = 3 

) 

DR 

< RE G* 2) 

ANp 

l R = 4 

) 

THEN 

1 1 : 

= 3 ; 

I F 

( ft E G = 1 ) 

and 

< 3x4 

) 

DR 

( R E G = 2 > 

A N p. 

< r> = 5 

) 

THEN 

1 1 : 

s 4 ; 

I F 

( R E G = 1 ) 

and 

<3 = 5 

) 

DR 

< R6G = 2) 

ANp 


) 

THEN 

1 1 : 

= s ; 

I F 

( REGsl ) 

and 

< 3X6 

) 

DR 

< ft E G = 2 ) 

A IN p 

< P = 7 

) 

THEN 

1 1 : 

= 6 ; 

1 F 

( R E G = 1 ) 

and 

<3=7 

) 

DR 

( R E G = 2 > 

AND 

< ft = 8 

> 

then 

1 1 : 

* 7 ; 

I F 

< R E G = 1 > 

and 

< 3 = 3 

) 

DR 

< ft E G= 2 > 

And 

(7 = 9 

) 

THEN 

1 1 : 

= 8 ; 

I F 

< R E G = 1 ) 

And 

< 3 = V 

) 

DR 

< R E G = 2 ) 

ANp 

3R=1 0) 

then 

1 1 : 

= 9 ; 

I F 

<REG=1 > 

A N D 

<3=10) 

DR 

< RFG = 2> 

ANp 

t ?e 1 1 

> 

then 

1 1 : 

= 10 ; 

I F 

( R€G = 1 ) 

And 

(3 = 11) 

DR 

( R fc G = 2 ) 

and 

l Pci 2) 

then 

I l t 

all; 

I F 

< R E G = 1 > 

and 

(3=12) 

dr 

< R £ G = 2 ) 

ANp 

C ? = 1 3 ) 

THEN 

1 1 : 

= 12 ; 


CASE II 0 F 


D 

I T E LN < ' 

ME * H 0 L E 

3 E ft I D D 

OF 

OBSERVATIONS' >; 

1 

W R I T E L \ < ' 

WINTERCFDR 

the 

HHOLt 

period 

OF 

OBSERVATIONS)*); 

? 

W ft I T E LN < 1 

S9RJMGCFDR 

the 

wHOlfc 

period 

OF 

OBSERVATIONS)'); 

5 

Wft I TE LN < ' 

$ J M M t « C F 0 a 

THE 

WHOLE 

period 

OF 

OBSERVATIONS)*),' 

!♦ 

W R I T E LN < ' 

A J T J vn < F OR 

the 

4 H 0 L t 

ptaioo 

OF 

OBSERVATIONS)'),' 

5 

W R I r E L N < ' 

4 I N T E a < 

F 

10.7 

= 

20o 

> ' ) ; 



8 

*i R I T E t \i < ' 

s a a i \ g c 

F 

10.7 

s 

200 

> ' > ; 



7 

WRIT E l N < ' 

S J ■* M F R C 

f 

10.7 

s 

2 Oo 

> ' > ; 



3 

w ft I T F L \ < ' 

A J T IJ M m ( 

F 

10.7 

s 

200 

) ' ) ; 



9 

w R I T E L N < ' 

4 I V T E a c 

F 

10,7 

= 

70 

> ' > ; 



1 D 

Wft I TE LN < ' 

S 9 a I N G c 

F 

10.7 

S 

70 

) ' ) ; 



1 1 

Wft I T E L\ < ' 

5 J M M F a C 

F 

10.7 

3 

70 

) ' > ; 



1 l 

wft I TE IS ( ' 

a j r u w n c 

F 

10.7 

s 

70 

> ' > ; 



END 

C * C A S E II 

* ) 









end; ( * NAME *) 

PROCEOJRE r A B L E ; 

C * RROCEOJRE PRINTS THF TABLE OF CALCULATED VALUES * ) 

(* OF DRIFT VELO-CIT/FS VX AND VY FOR EVERY *> 

C * HOUR FROM OOLT TILL 2 3 U FOR LATITUDES 5,10, IS, * ) 

<* . . . , 60 OEG. *) 

VAR I , J : integer; 

beg IN (* TABLE*) 

W R I T E L N < ' ZO'NAL VELOCITY VX 

WR I T E IN ; NR I T E ( ' '),-F0R J : * 1 TO 91 DO WRI TE C '**>,' WR I TE LN J 

WRITEC' FI L T » > I F 0 R J: = 1 TO 12 DO W « l T E C C J - 1 > : 7 ) ,' W R I T E L N »' 

WRITE(< ' > ; F 0 R J : a 1 TO 91 DO W R I T t < ' * * ) ; W R I T E L N J 

FOR J : =0 TO 10 00 

BEGIN 

J R I T E ( ' ',(60-J*5):2,' ' > ,' 

FOR I : = 0 TO 11 DO UR I T6 ( D X ( . 1 o- J , I , ) : ? : 1 ) 5 WR I TE LN 
end; 

WRITEC' ' ) J F 0 R J : = 1 TO 91 DO W R I T t ( * * * ) ; W R I T E L N ! 

wRiTEC' fi lT'>;for j : = i to 12 du wkitf<cj+iij:7);writfln; 

WRITEC' '),'FOR j:* 1 TO 91 DO vRlTfc(***);WRlTELN; 
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FDR J : = 0 TO ID DD 
BEGIN 

WRITER ' , ! SD - J *5 > 5 2 , ' '>? 

FDR I;sl2 TD ?J DO 4 R I T E < DX < , 1 0 - J , I , > J 7 J 1 ) ] WR I T E IN 
E V D / 

WRITEC » > ; F D R J * 3 1 TD 91 DD WRlTfe<'*f ) ;WRITEtNJWRITELN; 

U R I T E L N ! ' MERIDIONAL VELOCITT Vt < M / S > ' > J 

writeln; write!' p >;for J : * i to 91 do jrt t e ( ' * n ; wri te ln ; 

WRITE(* C I L r ' ) ; F 0 R J : = 1 TD 12 DO W^ITcC(J"1>*7)I w RT^ELNJ 

WRITEC* ' > J F D ft J I s 1 TD 91 DD WRlTt('*»>;WRlTELN; 

FDR J:=0 TO ID DD' 

BEGIN 

WRITE!' ',!5D-.J*5>:2 f ' ? > » 

FDR I : B o TD 11 OD wRlTE(0r(.lO-J,I.):7:1>IWRITELN 

e no ; 

WRITE!' ' > ; F D ft J : 3 1 TD 91 DD WRlTfc(»*f>;WRITELNJ 

WRITE!' FI L r ' > ; F 0 R J: = 1 TD 12 Do WKITFC(J*11) :7>fWRITELNJ 

WRITE!' ' ) ; F D R j: = 1 TD 91 DD W ft I T t ( ' * » ) J M ft I T E L N 

FDR J : = 0 TO ID D D 


BEGIN 

WRITE!' ',!0-J*5>:2 r ' ’ ) • 

FDR I : = 1 2 TD'23 DD W 9 I T E < D Y < . 1 0 - J i I 


) : 7: 1 > ; WRI TELN 


end; „ . „ , k 

WRITE!' 'If'FDR J : * 1 TD 91 DD WRlTfc(»*f);wrtlTELNjwRITELN 

end; ! * table *> 


BFGI* !* GLOBMDO *> 


v x : = D ; 'J 1 : = d ; 

for REG:*1 TO ? DD' 

BEGIN 

FDR I : S 1 TD 15 OD F D R J : = 1 TD f ft DD AR(.J,I.)l=OJ 
CASE REG OF 
1 : BEGIN 

!* INRJT the INITIAL data FOR E-REGIDN *> 
RESFT ! F A ft 1 > J 
FDR J : =1 TO ?2 DD 
BEGIN 

FDR I : 3 1 ro IS OD READ !FAR1,BFAR!.I,)>; 
FDR I ; 3 1 TD 15 DD AR<.J,I.>: = 3FAR(.I f > 
END 


? : 


end; 

BEGIN 

! * 


! * R E G1 * ) 

INRJT THF INITIAL DATA FOR F ” R E G 1 D N 


* > 


reset !Far?>; 

FOR J : * 1 TO ? B DD 
BEGIN 

FDR I J s 1 TD 15 DD REAp !FAR2,BFAR(.I,)>J 
FDR I ; = 1 TO 15 DD A K < . J , I . > : s B F A R < , I , > 


END 


END!* R E 32 * > 

end; ( * CASE * > 


CASE REG OF 

1 : P 1 AX ; = 1 2 ? 

2 : P 4 A X ; * 13 

e^o; 
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FOR P:s1 TO MM 00 
0 E G I * 

P A 3 E ; 

NA'IECREG,:*); 

F I : = 5 ; 

for J:=o ro io do 

3EGI M 

f I ! * F I *■ 5 ; 

FOR I : =0 TO 23 00 
3E3I * 

L T l = I ; 

ORIRR(PEG» & fFIfLT,Vx,vT); 

E V 0 <* 1 *) 

6 n o ; ( * J * ) 

TA UE 

END ( * P * ) 

F N 0 < * REG * ) 

END. C * GLOB ^00 * > 
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